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•  I 


Dlit 


For  the  time  evolution  of  a  Hele-Shaw  interface  described  by  a  conformal  map  z{C,t) 
that  maps  a  unit  circle  (or  a  semi-circle)  in  the  C  plane  into  the  viscous  fluid  flow  region 
in  the  physical  2-plaiie,  we  present  results  on  the  motion  of  singularities  outside  the  unit 
circle. 

For  zero  surface  tension,  we  extend  earlier  results  to  show  that  for  any  initial  condition, 
each  singularity  of  present  initially  in  |C|  >  1  continually  approaches  the  interfacial 

boundary  |C|  =  1  without  any  change  of  form.  However,  depending  on  the  singularity 
type,  it  may  or  may  not  impinge  |C|  =  1  hi  finite  time.  Under  some  assumptions,  we  give 
analytical  evidence  to  suggest  that  the  ill-posed  problem  in  the  physical  domain  |C|  <  1 
can  be  imbedded  in  a  well-posed  problem  in  |C1  >  1.  We  present  a  numerical  scheme  to 
calculate  such  solutions. 

For  each  initial  singularity  of  certain  type,  which  in  the  absence  of  surface  tension 
would  have  merely  moved  to  a  new  location  Cs(^)  at  time  t  from  an  initial  Cs(0)>  we  find 
an  immediate  transformation  of  the  singularity  structure  for  nonzero  surface  tension  H; 
however,  for  0  <  B  <<  1,  surface  tension  effects  on  this  singularity  are  limited  to  a 
small  ’inner’  neighborhood  of  Cs(0  when  t  <<  i.  Outside  the  inner-region  but  for 
1^  —  <<  1,  the  singular  behavior  of  the  zero  surface  tension  solution  still  persists 

for  z{(^,t).  On  the  other  hand,  for  each  initial  zero  of  z^,  which  for  surface  tension  B  =  0 
remains  a  zero  of  at  a  location  Co(0  different  from  Co(0)>  surface  tension  effects  spawns 
new  singularities  that  move  away  from  Co(t)  approach  the  physical  domain  |C|  =  1. 
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1.  Introduction 

The  problem  of  a  less  viscous  fluid  displacing  a  more  viscous  fluid  in  a  Hele-Shaw 
cell  has  received  considerable  attention  over  the  last  few  years  since  the  original  work  by 
Saifman  &  Taylor  (1958).  Reviews  by  Saffman  (1986),  Bensimon  et  al  (1986),  Homsy 
(1986)  svunmarize  the  state  of  affairs  as  of  flve  years  back.  Since  then,  there  has  been 
considerable  work  that  have  been  reviewed  from  different  perspectives  by  Pelce  (1988), 
Kessler  &  Levine  (1988),  Howison  (1991)  and  Tanveer  (1991).  This  intense  activity  related 
to  a  relatively  simple  experiment  has  been  motivated  by  the  mathematical  analogies  of  the 
Hele-Shaw  to  viscous  displacement  in  a  porous  medium,  growth  of  a  needle  crystal  in  an 
undercooled  melt  and  the  morphological  instability  and  resulting  features  in  directional 
solidiflcation.  Indeed,  many  of  the  generic  results  for  these  apparently  diverse  areas  were 
first  discovered  in  a  Hele-Shaw  flow.  In  addition  to  the  relative  simplicity  of  the  equations, 
the  flow  in  this  geometry  is  accessible  to  experiments  making  it  possible  to  compare  theory 
with  experiment. 

Most  of  the  work  to  date  has  been  on  steady  states  and  their  linear  stability  (see  recent 
reviews  cited  above  for  extensive  bibliography).  In  an  actual  experiment  (Saffman  &  Taylor 
(1958),  Tabeling  et  al  (1986),  Maxworthy  (1987),  Arneodo  et  al  (1989)  the  interface  does 
not  evolve  into  a  steady  state  when  surface  tension  effects  are  very  small.  Indeed,  through 
a  sequence  of  instabilities,  a  highly  convoluted  interfacial  pattern  forms  that  appears  to 
be  fractal  (Arneodo  et  al  (1989))  over  some  range  of  length  scales.  Further,  the  precise 
interfacial  shape  appears  to  be  extremely  sensitive  to  initial  conditions  though  the  overall 
statistics  of  the  observed  pattern  is  not.  The  theoretical  understanding  of  these  features 
is  lacking  to  a  great  degree.  This  paper,  we  hope,  is  an  important  step  in  that  direction. 

The  time  evolution  problem  is  studied  in  this  paper  is  for  a  highly  idealized  boundary 
condition,  originally  developed  by  by  Mclean  &  Saffman  (1980)  (we  call  it  the  MS  bound¬ 
ary  conditions),  where  important  three  dimensional  effects  observed  in  the  experiment 
(Tabeling  et  al  (1986))  are  ignored.  For  steady  flow,  realistic  two  dimensional  boundary 
conditions  incorporating  three  dimensional  thin  film  effects  were  discussed  by  Saffman 
(1982)  and  later  derived  by  Park  &  Homsy  (1985)  and  Reinelt  (1987)  (we  call  it  the  SPHR 
boimdary  conditions).  Numerical  (Schwartz  &  DeGregoria  (1987),  Reinelt  (1987),  Sarkar 
&  Jasnow  (1987))  and  analytical  calculations  (Tanveer,  1990)  show  important  quantitative 
differences  in  solutions  corresponding  to  the  MS  and  SPHR  boundary  conditions.  Over 
some  range  of  experimental  control  parameters,  qualitative  differences  also  exist  though  in 
other  cases  the  qualitative  conclusions  zu:e  about  the  same.  Further,  analytical  calculations 
on  steady  state  finger  (Tanveer,  1990)  suggest  that  the  mathematical  structure,  in  a  broad 
sense,  are  similar  in  these  two  cases  despite  the  relative  complexity  of  the  SPHR  conditions. 
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Here,  we  proceed  with  the  belief  that  for  the  time  evolving  flow,  the  solutions  correspond¬ 
ing  to  the  much  simpler  MS  boundary  conditions  are  relevant,  atleast  qualitatively,  over 
some  range  of  control  parameters,  as  is  true  for  the  steady  state. 

For  the  rectilinear  geometry,  we  take  z{(i,i)  to  be  the  conformal  map  that  maps  the 
semi-circle  in  Fig.  1  into  the  physical  flow  domain  (Fig.  2)  in  the  z  plane  such  that  ^  =  0 
corresponds  to  2;  =  -|-oo  and  ^  =  ±1  correspond  to  points  A  and  B  (see  Fig.  2),  where  the 
interface  meets  the  channel  side  walls.  It  is  clear  that  one  can  decompose 

^'(C,  t)  =  --  /nC  +  i  +  /(C,  *)  (1-1) 

K 

where  /  is  analytic  in  the  unit  semi-circle  with 

Im  /  =  0  (1-2) 

on  the  real  diameter  of  the  unit  circle  corresponding  to  the  geometric  condition  that  the 
walls  of  the  channel  correspond  to  the  real  diameter  of  the  unit  semi-circle.  We  will  assume 
/  is  continuous  up  to  the  real  diameter  including  (^  =  0  .  Further,  we  will  assume  that 
the  shape  of  the  extended  Anger  obtained  by  reflection  about  each  of  the  two  side  walls  is 
smooth  implying  that  z  and  hence  /  is  analytic  and  Z(^  ^  0  on  the  semi-circular  arc 
I  Cl  =  !•  iFrom  Schwarz  reflection  principle,  /  is  analytic  and  ^  0  for  |C|  <  1  . 
We  decompose  the  complex  velocity  potential: 

W(C,  i)  =  --  In  C  +  i  +  w(C,  t)  (1.3) 

TT 

where  the  fluid  velocity  at  infinity  is  assumed  to  be  unity  without  any  loss  of  generality 
(equivalent  to  appropriate  nondimensionalization).  It  is  clear  that  the  condition  of  no  flow 
through  the  walls  imply 

Im  u;  =  0  (1-4) 

on  the  real  diameter  (-1,  1)  of  the  unit  semi-circle.  As  is  physically  reasonable,  a;  will  be 
Eissumed  to  be  continuous  up  to  the  real  diameter  including  C  =  0  .  Further,  it  will  be 
assumed  that  u  is  analytic  on  the  semicircle  corresponding  to  the  assumption  of  smooth 
flow  at  the  interface.  ^From  Schwarz  reflection  principle,  (1.4)  implies  that  u)  is  analytic 
in  Id  <  1  • 

Second,  for  the  radial  geometry,  following  Howison  (1985,1986ab),  the  conformal  map 
2(C,  t)  that  maps  the  interior  of  the  unit  circle  in  the  C  plane  into  the  physical  flow 
domain  (Fig.  1)  such  that  ^  =  0  is  mapped  to  oo  can  be  decomposed  as 

KC,  0  =  ^  +  MC.  t)  (1.5) 
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where  the  remaining  degree  of  freedom  in  the  Riemann  mapping  theorem  is  used  by  re¬ 
quiring  a{t)  to  be  real  and  positive.  We  will  assume  that  the  ili•,el^ace  is  smooth  and  so 
k{^,  t)  is  analytic  and  Z(^  ^  0  everywhere  in  |C|  <  1-  c.m  be  expressed  as  as 

power  series  convergent  for  |C|  <  1  : 

OO 

KM  =  Y,  *"(')  <"  (1-6) 

n=0 

Further,  the  complex  velocity  potential  t)  can  be  decomposed  as: 

W((,  t)  =  -  -In  (  +  u,((,  t)  (1.7) 

TT 

where  the  injection  rate  Q  =  4  without  any  loss  of  generality.  This  choice  is  made  for 
convenience  of  presentation  so  that  (1.3)  and  (1.7)  have  the  same  form.  For  smooth  flow 
at  the  interface,  t)  can  be  assumed  to  be  analytic  in  |^|  <  1  and  has  a  convergent 
series  representation  in  that  domain: 

OO 

n=0 

Note  that  unlike  the  channel  geometry,  w  and  k  do  not  have  vanishing  imaginary  part 
on  the  real  axis  within  the  unit  circle  except  if  they  are  so  initially  (corresponds  to  certain 
conditions  of  symmetry  of  initial  conditions). 

In  either  Hele-Shaw  cell  geometries,  the  boxmdary  condition  (Mclean-Saffman,  1980) 
on  the  finger  boundary  that  the  difference  of  pressure  on  the  two  sides  is  balanced  by 
surface  tension  times  curvature  implies  that  on  |C|  =  1 

Reuj  =  — Re  I  +  (1-9) 

Pel 

where  for  the  channel  geometry,  B  =  i2fiVa^  ’  ^  is  the  surface  tension,  h  is  the 

gap  width  and  2a  is  the  cell  width,  the  viscosity  of  the  more  viscous  fluid,  V  the 
velocity  of  fluid  at  oo  .  For  the  radi2il  Hele-Shaw  geometry,  we  define  B  =  ,  where 

Q  is  the  injection  rate  and  tt  is  the  initial  area  of  the  blob  of  less  viscous  fluid.  The 
viscosity  of  the  less  viscous  fluid  and  the  thin  film  effects  have  been  neglected  here  for  the 
sahe  of  simplicity.  The  kinematic  boundary  condition  that  corresponds  to  no  fluid  flow 
through  the  interface  (see  Saffman  (1959)  or  Richardson  (1972)  for  details)  implies 


(1.10) 
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on  I  Cl  =  1  .  Note  that  for  the  channel  geometry,  (1.9)  and  (1.10)  hold  on  the  lower  half 
unit  semi-circiilax  boimdary  as  well.  The  mathematical  problem,  therefore,  is  to  determine 
analytic  functions  u>  and  /  (or  k  for  the  radial  geometry)  and  therefore  W{C,t)  and  z{C,t) 
satisfying  (1.9)  and  (1.10)  on  the  circular  arc,  when  ^(C,0)  and  W(C,0)  axe  appropriately 
specified. 

Almost  all  the  work  to  date  on  this  initial  value  problem*  belong  to  either  of  two 
categories:  numerical  computations  for  B  not  too  small  or  analytical  work  for  B  =  0. 

For  B  0,  the  numerical  computations  by  Trygvasson  &  Aref  (1983),  DeGregoria  & 
Schwartz  (1986),  Bensimon  (1986)  and  Meiburg  &  Homsy  (1988)  have  shed  some  light  on 
the  process  of  finger  competition  that  eventually  leads  to  a  single  steadily  moving  fingf^r. 
However,  when  surface  tension  parameter  is  smaller  than  some  critical  value,  the  finger  is 
itself  unstable.  The  calculations  of  linear  stability  by  Kessler  &  Levine  (1985,  1986)  and 
Tanveer  (1987b)  suggests  that  there  exists  a  linearly  stable  branch  for  any  nonzero  surface 
tension.  Based  on  numerical  experimentation,  Bensimon  (1986)  concluded  that  observed 
instability  in  the  simulation  is  a  nonlinear  mechanism  where  the  threshold  amplitude  for 
destabilization  decreases  with  surface  tension.  This  is  plausible,  as  a  discrete  set  of  steady 
fingers  (VandenBroeck,  1983)  coalesce  as  B  0  and  all  but  one  of  these  branches  are 
unstable  (Kessler  &  Levine  (1987),  Tanveer  (1987b),  Bensimon  et  al  (1987)).  Despite  the 
understanding  on  the  onset  of  various  instabilities,  a  full  understanding  of  the  complicated 
time  dependent  pattern  beyond  the  onset  stage  is  difficult  using  conventional  simulation 
with  boundary  integral  methods  since  widely  disparate  scales,  both  in  space  and  time, 
need  to  be  resolved  accurately  as  B  is  made  progressively  smaller. 

For  B  =  0,  it  is  clear  that  the  initial  value  problem  posed  in  (1.9)  and  (1.10)  sim¬ 
plify  considerably  since  u;  =  0  in  this  case  and  therefore  W{C,t)  is  known  (from  (1.7)). 
Gustaifson  (1984,  1987)  has  rigorously  prove  a  the  existence  of  solution  for  finite  time  for 
analytic  initial  data  on  the  real  axis  (the  physical  domain  in  his  formulation).  Earlier, 
Galin  (1945)  and  Polubarinova-Kochina  (1945)  considered  the  mathematically  identical 
problem  of  Darcy  model  for  ground  water  flows  and  devised  analytical  techniques  for  ob¬ 
taining  exact  solutions.  These  were  apparently  well  known  in  the  Russian  literature  (see 
Howison  (1991)  &  Hohlov  (1990)).  Exact  solutions  due  to  Saffman  (1959),  Howison  (1985, 
1986ab)  and  Shraiman  &  Bensimon  (1985)  can  be  seen  as  applications  of  these  techniques 
though  results  appear  to  have  been  obtained  without  knowledge  of  the  Russian  work. 
Howison  (1990)  summarizes  the  relation  between  the  different  techniques.  The  SaJFman- 
Howison- Shraiman- Bensimon  (SHSB)  exact  solutions  are  limited  to  initial  conditions  for 
which  ^^(C>  0)  is  either  a  polynomial  or  a  superposition  of  poles.  In  these  cases,  the  number 

*  Other  authors  have  used  differing  but  equivalent  formulation 


of  poles  is  known  to  be  preserved,  as  is  the  number  of  zeros;  only  their  locations  change. 
Within  this  class  of  solutions,  it  is  known  that  poles  never  hit  the  physical  domain  bound¬ 
ary  1<^|  =  1  in  finite  time  but  zeros  may  or  may  not.  In  the  case  when  a  zero  hits  |^|  =  1, 
a  zero  angled  cusp  protruding  into  the  viscous  fiuid  forms  at  the  physical  interface.  The 
mathematical  solution  does  not  make  physical  sense  beyond  that  time.  Howison  (1986a) 
was  further  able  to  show  that  that  one  can  start  with  two  different  initial  conditions  that 
axe  arbitrarily  close  in  the  physical  domain  |Cl  <  1  (in  any  Sobolev  norm),  yet  after  a 
finite  time,  the  physical  boundary  develops  a  cusp  in  one  case  while  it  remains  regular  in 
the  other.  In  this  sense,  it  is  known  that  the  zero  surface  tension  evolution  problem  is  ill 
posed.  Richardson  (1972)  showed  certain  invariance  properties  of  the  Schwartz  function 
of  a  domain  for  the  related  problem  of  a  more  viscous  fiuid  displacing  a  less  viscous  one. 
In  our  formulation,  Richardson’s  result  implies  that  singularities  of  ^(Cj*)  in  the  finite  ( 
plane  other  than  that  at  C  =  0i  arise  due  motion  of  an  initial  singularity  of  the  same  form 
(i.e.  for  instance,  a  branch  point  of  one  half  power  remains  so  at  later  times)  of  2((^,  0)  in 
|C|  >  1.  However,  no  explicit  solution  are  known  if  Z((C,  0)  contains  branch  point.  Richard¬ 
son’s  suggestion  about  approximating  them  by  a  polynomial  is  invalid  in  our  case  as  we 
shsdl  see  that  the  later  time  evolution  depend  critically  on  the  initial  singularity  structure 
in  j^l  >  1.  The  only  known  closed  form  solutions  with  branch  points  are  the  self-similar 
solutions  in  time  for  a  radial  or  wedge  geometry  recently  foimd  by  Thome  et  al  (1988), 
Ben  Amar  (1991)  and  Tu  (1991).  Further,  due  to  the  ill-posedness  of  the  problem  in  the 
physical  domain  |^|  <  1,  conventional  numerical  simulation  using  boundary  integral  meth¬ 
ods  for  instance  is  not  possible  due  to  vmcontrolled  growth  of  roxmd  off  errors,  unless  one 
uses  some  artificial  dissipation  or  a  filtering  procedure  (see  Krasny  (1986)  for  applications 
to  another  ill-posed  interfacial  problem).  Such  procedures  appear  to  be  computationally 
expensive  and  may  not  be  practical  when  the  interface  becomes  highly  convoluted  as  is 
the  case  in  the  interesting  stages  of  Hele-Shaw  evolution. 

For  nonzero  B,  it  is  expected  that  the  initial  value  problem  in  the  physical  domain 
will  be  well-posed.  However,  the  problem  is  considerably  more  intractable  analytically, 
because  of  complications  due  to  the  unwieldly  curvature  term  on  the  right  hand  side  of 
(1.9).  As  stated  earlier,  numerical  computations  in  the  physical  domain  become  rather 
difl&cult  in  the  dynamically  interesting  case  of  small  B.  Duchon  &  Robert  (1984),  while 
working  on  a  mathematically  equivalent  problem,  have  proved  the  existence  of  solutions 
for  short  times.  However,  it  is  not  known  if  solutions  exists  for  all  times  for  any  nonzero 
B.  Physically,  one  might  expect  this  to  be  the  case,  though  this  remains  to  be  shown. 
Studying  the  problem  perturbatively  about  any  known  solution  at  H  =  0  is  fraught  with 
difficulties.  If  So{t)  denotes  the  solution  operator  that  maps  the  initial  condition  2(C>0)  in 
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ICl  <  1  into  the  6  =  0  solution  in  |^|  <  1,  ill  posedness  immediately  implies  that  the 

operator  So(t)  is  unbounded.  For  this  reason,  there  is  no  reason  to  expect  under  a  small 
perturbation  caused  by  a  small  nonzero  6,  the  corresponding  solution  operator  S(t)  will 
be  close  to  6o(t).  Any  physical  argument  to  the  contrary  is  unacceptable  as  a  physical  flow 
always  corresponds  to  nonzero  B  and  one  cannot  use  physical  argument  to  extrapolate  the 
6  =  0  result.  Despite  this,  in  an  effort  to  make  some  headway,  analytical  efforts  have  been 
made  by  Lacey  et  al  (1988)  by  assuming  that  for  small  surface  tension,  surface  tension 
effects  are  important  only  when  the  interface  curvature,  as  predicted  by  z°,  is  large.  As 
the  authors  themselves  point  out,  this  assumes  that  the  actual  interface  curvature  is  not 
very  different  from  what  is  predicted  by  2°,  the  zero  surface  tension  solution.  Recently, 
Dai  et  al  (1991)  carry  out  a  numerical  simulation  of  the  initial  value  problem  in  the 
physical  domain  for  the  small  nonzero  6  in  the  radial  geometry  when  Z(  is  either  a  cubic 
polynomial  or  a  combination  of  three  poles  with  a  three  fold  rotational  symmetry  and 
deduce  the  nature  of  singularities  in  the  unphysical  domain  using  a  Domb- Sykes  plot.  For 
such  initial  conditions,  exact  solutions  are  known  within  the  Saffman-Howison-Bensimon- 
Shraiman  family  when  6  =  0.  Comparison  of  Dai  et  al  (1991)  calculations  with  these 
suggest  that  initial  poles  of  Z(^  remain  unaffected.  With  polynomial  initial  conditions,  it 
is  found  that  an  initial  zero  splits  into  what  they  term  as  ’two  pole-like  singularities’.  In 
this  case,  it  is  found  that  difference  with  the  exact  6  =  0  solution  becomes  significant  even 
when  the  interface  curvature  based  on  2°(^,<)  is  <<  This  appears  to  suggest  that  the 
Lacey  et  al  (1988)  hypothesis  may  not  be  correct.  Recently,  Constantin  &  Kadanoff  (1990, 
1991)  have  derived  local  partial  differential  equations  for  in  the  region  |^|  >  1  for 

the  radial  geometry  by  analytically  continuing  (1.9)  and  (1.10)  for  the  radial  geometry, 
assuming  that  the  singularities  are  very  far  off  from  the  physical  plane  so  that  the  physical 
interface  is  nearly  circular.  Their  localization  procedure  is  similar  in  spirit  to  the  equations 
derived  by  Caflisch,  Orellana  &  Siegel  (1990)  for  the  zero  surface  tension  case.  Constantin 
&  Kadanoff  (1991)  and  Howison  (1991)  have  also  found  special  similarity  solutions  that 
show  the  immediate  transformation  of  an  initial  singularity  or  a  zero  of  2^  into  possible 
-4/3  singularities;  however  their  relevance  to  the  overall  physical  dynamics  has  not  been 
addressed.  Constantin  &  Kadanoff  (1990)  also  prove  the  existence  of  solutions  to  these 
localized  equations  in  an  unphysical  domain  that  shrinks  in  time  like  .  Since  the 

physical  domain  is  |^|  <  1  (in  their  formulation,  |^|  >  1  ),  their  results  do  not 

guarantee  existence  of  solutions  for  all  times  even  with  surface  tension.  However,  their 
equations  become  invalid  when  the  actual  interfacial  shape  deviates  significantly  from  a 
circle,  long  before  a  singularity  actually  can  come  close  to  |C|  =  1-  Thus  the  behavior  of 
their  solution  need  not  reflect  the  behavior  actual  solutions  that  satisfy  (1.9)  and  (1.10). 
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We  like  to  iinderstaiid  the  Hele-Shaw  dynamics  for  small  nonzero  B  perturbatively  by 
exploiting  the  simplicity  of  the  equations  when  B  =  0.  There  are  two  major  hurdles  in 
accomplishing  this.  First,  the  information  on  the  B  =  0  problem  is  not  complete.  We  do 
not  have  exact  solutions  for  a  general  initial  condition.  Nor  is  it  known  how  to  compute 
such  solutions  effectively.  Second,  B  =  0  problem  is  ill  posed  in  the  physical  domain  |C|  <  1, 
as  demonstrated  by  Howison  (1985,  1986ab).  Thus,  as  discussed  before,  one  cannot  have 
any  confidence  that  a  zero  surface  tension  solution  ^“(C,  t)  in  the  physical  domain  ICl  <  1 
is  necessarily  relevant  as  B  — »  0.  In  this  paper,  we  address  these  two  problems. 

First,  we  show  that  the  analytically  continued  equation  for  B  —  0  has  the  feature  that 
all  information  in  |C|  >  1,  flows  inwards  towards  the  interface  boundary  |C|  =  1-  In  partic¬ 
ular,  for  arbitrary  initial  conditions,  singularities  present  in  |C|  >  1  must  always  approach 
the  physical  domain  boundary  |^|  =  1,  though  it  need  not  actually  impinge  it  in  finite 
time.  Assuming  that  that  the  speed  of  a  zero  of  Z((C,  f)  is  bounded  for  certain  classes  of 
initial  conditions,  we  present  analytical  arguments  to  suggest  that  the  initial  value  prob¬ 
lem  in  Id  '>  1  obtained  by  solving  the  analytically  continued  equations  with  analytically 
continued  initial  condition  is  well  posed  in  the  unphysical  domain  for  a  finite  time.  We 
also  argue  that  the  solution  obtained  in  this  way  for  <  >  0  is  the  analytical  continuation 
of  the  physical  solution  in  |d  <  1  across  the  unit  circle.  Using  these  results,  we  suggest  a 
numerical  method  to  find  solutions  for  arbitrary  initial  conditions,  including  ones  where 
has  a  set  of  branch  points  of  specific  types.  Because  of  the  well-posedness  of  the  underly¬ 
ing  formulation  in  |d  >  1,  the  computation  can  be  performed  without  growth  of  round  off 
errors,  a  problem  that  plagues  computation  of  solutions  to  ill-posed  problems.  Recently, 
G.R.  Baker  (private  communication)  has  successfully  implemented  this  scheme.  We  then 
present  analytical  evidence  that  a  singularity  where  z^  -Bo(<)(C  -  G(0)  does 
not  impinge  the  physical  domain  when  0  When  0  <  /?  <  |,  our  evidence  suggests 

that  the  singularity  impinges  |d  =  1  in  finite  time.  Under  some  assumptions  that  we  are 
unable  to  check,  the  same  appears  to  be  true  for  <  0.  Our  results  also  suggests  that  the 
ill-posedness  in  the  physical  domeun  for  B  =  0  is  actually  in  the  sense  of  Hadamard. 

Having  found  a  well  posed  formulation  in  the  domain  |C|  ^  1  allows  us  to  address  the 
second  serious  hurdle  mentioned  above.  In  this  extended  domain,  we  can  expect  that  the 
solution  for  nonzero  is  close  in  some  sense  to  z°(C,t)  and  therefore  can  be  studied 

perturbatively.  Indeed,  a  minor  modification  of  the  numerical  method  suggested  for  solving 
the  zero  surfzice  tension  solution  z°(C,f)  can,  in  principle,  be  used  to  find  each  term  of  the 
regular  perturbation  expansion  in  |C1  >  1: 

ziM  =  -h  Bz*(C,0  +  ...  (1.11) 

Examination  of  the  analytically  continued  equation  m  |(|  >  1  shows  immediately  that  the 
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regular  perturbation  expansion  expansion  will  break  down  near  each  point  Co(t)  where  = 
0  amd  singularities  that  include  branch  points  CsCO  which  2°  ~  Bo(t)(C  —  Ca(0)  ^ 
with  0  <  2.  Even  for  B  independent  initial  conditions,  the  perturbation  expansion  (1.11) 

^”4  /3  ^ 

becomes  invalid  at  any  singular  point  (^p  where  z,^  ^(f)[c  -  Cp(<)]  .  with  A(0)  =  0 

and  A{t)  scaling  as  some  positive  power  of  B.  As  we  shall  see  later,  such  singular  points 
are  created  at  the  initial  instant  of  time  by  surface  tension  effects.  When,  the  perturbation 
expansion  (1.1)  breaks  down,  it  becomes  necessary  to  consider  the  inner  equation  obtained 
by  appropriately  rescaling  the  dependent  and  independent  variables. 

With  certain  assumptions,  we  find  that  for  initial  conditions  that  are  independent  of 
B,  appropriate  inner-outer  matching  can  be  carried  out  for  0  <  0  <  2.  We  find  that  for  a 
branch  point  or  a  simple  pole  singularity  of  2,j((^,0)  where  2(^(C,0)  ~  jBo(O)  (C  —  Ca(0))  > 

for  a  nonzero  B,  the  singularity  is  immediately  transformed.  However,  for  0  <  B  << 
1  and  t  <<  -g,  the  location  of  the  newly  created  singularities  is  within  a  small  inner 
neighborhood  of  (s(t)  (singularity  of  z^),  where  |C  —  C»(0I  =  0(Ba(*-^> ).  The  singular 
behavior  ((  —  0(0)”^  reflected  in  the  actual  solution  z^  when  1  >>  |C  —  Ca(0l 

BSTj^T.  Thus,  in  the  outer  asymptotic  sense,  the  solution  z{(^,t)  behaves  like  2®((,t)  near 
Ca{t).  When  |Cs(f)|  — »  1,  as  must  happen  for  every  singularity  at  the  late  stages,  the  physical 
interface  can  become  quite  distorted  locally.  In  this  final  stage,  our  analysis  is  restricted 
to  5  <  /?  <  1.  We  find  that  the  governing  equations  when  |Cs(<)l  ~  1  is 

nonlocal  in  an  inner  scale  where  ^  —  C4(0  “  We  have  not  attempted  to 

solve  this  inner  equation;  however,  the  in  this  inner  neighborhood  scales  as 
which  is  large  for  small  B.  We  conclude  that  when  a  singularity  of  this  kind  is  present  in 
the  initial  data  in  (^|  >  1,  it  will  eventually  cause  locally  large  localized  distortions  of  the 
interface  when  lCs(0l  1- 

For  an  initial  zero  of  2<;(C,0)  that  is  independent  of  B,  we  find  that  while  2°(^,t) 
would  have  merely  moved  the  zero  to  a  new  location  Co(0)  effect  of  surface  tension 
is  create  new  singularities  immediately.  The  possibility  of  these  singularities,  where  z^ 
has  a  -4/3  power  singularity  just  like  the  ones  for  an  initial  branch  point  singularity  was 
recognized  by  Constantin  &  Kadanoff  (1991)  and  Howison  (1991).  We  actually  confirm 
their  existence  and  calculate  their  locations  at  early  times.  Further,  unlike  the  case  of  initial 
branch  points,  where  such  singularity  locations  are  within  an  inner  region  around  C*(0i 
which  masks  the  effects  of  the  -4/3  singularities,  in  this  case,  each  of  the  -4/3  singularities 
move  away  from  Co(0-  ^  =  ^(1)»  these  singularities  are  clustered  over  a  B^^® 

scale.  However,  the  effect  of  each  of  these  singularities  Cp(0  ^  small  B^^^* 

neighborhood  around  each  of  these  points  beyond  which  the  solution  matches  to  2°(Cf). 
In  this  region  of  singularity  cluster  created  by  an  initial  zero  Co(0),  there  is  also  a  point 
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Cji(<)>  where  the  perturbation  expansion  (1.11)  breaks  down  because  is  singular  and 
higher  order  terms  2^,  2^,  etc.  are  progressively  more  singular.  However,  this  is  not  a  true 
siij^aiarity  of  2^,  as  the  solution  is  smoothed  out  in  ^m  inner  scale  of  extent  where 

the  deviation  of  Z(^  from  2°  is  at  best  This  inner  regions  around  each  Cp(0 

C*i(0  do  not  overlap.  Our  results  here  become  invcdid  when  any  of  |Cp(01  |Cji(OI  very 
close  to  1,  as  in  that  case  new  scalings  of  the  inner  variables  become  necessary  that  lead 
to  new  nonlocal  equations  that  are  yet  to  be  einalyzed.  We  are  of  the  opinion  that  Dai 
et  al’s  (1991)  direct  numerical  computation  reflecting  the  creation  of  new  singularities  out 
of  an  initial  zero  is  reflected  in  our  findings,  though  the  detected  ’pole-like’  singularity  in 
the  numerical  calculations  may  be  an  artifact  of  the  clustering  of  the  singular  ities  over  a 
52/9  scale.  One  of  the  immediate  implications  of  our  findings  is  that  it  is  possible  for  2^ 
to  deviate  from  2°  by  0(1)  or  larger  even  at  points  where  2®  is  neither  zero  nor  singular. 

Our  results  for  small  nonzero  B  suggests  that  we  can  start  with  smooth  physical 
interface  initally;  at  later  times  the  interface  can  get  highly  convoluted  over  a  local  length 
scale,  as  the  effect  of  singularities  either  present  initially,  or  those  created  from  an  initial 
zero  of  2^  continually  move  towards  |^|  =  1.  Since  the  location  of  initial  singularity  and 
zeros  is  highly  sensitive  to  arbitrarily  small  changes  of  z{(^,  0)  in  the  physical  domain  jCl  <  1, 
there  will  be  highly  sensitive  dependence  of  observed  features  on  initial  conditions. 

2.  Analytically  continued  equations  in  |C|  >  1 

Using  Poisson  integral  formulae  that  relates  a  harmonic  function  and  its  conjugate  to 
its  boundary  value  on  the  unit  circle,  from  (1.9)  for  |C|  <  1  , 

u  =  -B  h  (2.1) 


where 


h  =  ^  i  rfC'  p 
27^1  7|ci=i  LC' 


+  c 


n 


HO 


where 


2c1/2(C',  t) 


I2  = 


1  ^CC  (C',0  ,  1 


1  +  -  c 

2  2c(C,f) 


+  2 


t) 

1  %  (J,  0 

C  kCc 


I,  =  - 


(2.2) 


(2.3) 


(2.4) 


where  for  the  channel  geometry. 


2(C,0  =  2(C,0 


(2.5) 
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and  for  the  radial  geometry: 


i(C,  <)  =  ^  +  E  KW  < 

^  n=0 


(2.6) 


where  superscript  *  denotes  complex  conjugate.  By  a  standard  method  of  analytic  con¬ 
tinuation  through  contour  deformation  (see  Tanveer,  1987a  for  an  example),  for  |^(  >  1  , 


U} 


-Bh 


B  h 


1 


2  +  C  —  +  /3 


(2.7) 


Alternatively,  (2.7)  follows  from  (2.1)  on  using  a  variation  of  the  well  known  Plemelj 
formula  (see  for  example  Carrier,  Krook  &  Person  (1983)).  Note  the  choice  of  branch  of 
half  power  appearing  in  (2.3)  and  (2.7)  should  be  consistent  with  |  =  on 

c  =  e-  . 

The  analytical  continuation  leading  to  (2.7)  can  alternately  be  done  by  noting  that 
on  ^  =  e’"  ,  an  analytic  function  G(C)  in  |C|  <  1  with  representation: 


0(0  =  E  9”  (2.8) 

n=0 

satisfies  the  relation:  2  Re  G  =  G{(,t)  -|-  G(^,  t)  ,  where 

G(c.<)  =  E  s;  c”  (2.9) 

n^O 

Applying  this  to  (2.9),  we  obtain  (2.7)  with 

1  ^'(i,  <)  (2.10) 

This  provides  an  alternate  expression  for  Ii  besides  (2.2). 

Using  Poisson’s  integral  formulae  for  the  harmonic  function  Re  [^]  ,  with  appropri¬ 
ate  choice  of  imaginary  constant  for  |^|  <  1  ,  (1.10)  implies: 


Zt  =  C  z<;  U 


(2.11) 


where 

h  (C,  0 


47ri 


[C  +  C'l 

C' 

Ic'  -  cj 

CWdC,  0  . „  ,  /6«',  0 

\{C',  0  +  7:77rT^ 


^c(C',  OJ 


(2.12) 
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where 


II 

(2.13) 

,  1  0 

'  C  z<(^  t) 

(2.14) 

By  urmg  analytic  continuation  through  contoxir  deformation,  for  |C| 

>  1  we  obtain: 

^4  +  C"  /5  +  C  le 

(2.15) 

Alternately,  if  we  use  the  procedure  of  the  last  paragraph,  we  can  derive  (2.15)  provided 

II 

(2.16) 

Equation  (2.16)  provides  an  alternate  expression  for  I4  (besides  (2.12))  for  |(  |  >  1  . 

Combining  (2.7)  and  (2.15),  we  obtain  one  nonlinear  integro-differential  equation  for 

z  for  1^1  >  1  : 

r,  =  ,.  r,  +  „  +  8,3  +  -  1 

+  8^  (2.17) 

where 

II 

(2.18) 

92  =  Ch  +  Ch 

(2.19) 

93  =  — h 

(2.20) 

94  =  —h<i  h  (2  +  h)  —  h  h  C  h(^ 

(2.21) 

95  =  I5  +  2  -^2  h  h 

(2.22) 

97  =  — I2  h 

(2.23) 

It  may  be  noted  that  each  of  q\  through  qj  are  emalytic  functions  of  (  outside  the  vmit 
circle  since  they  involve  ^  and  the  evaluated  on  or  inside  the  unit  circle  where 
they  are  analytic. 

3.  The  case  of  zero  surface  tension 

It  is  clear  from  that  when  B  =  0  ,  uf  =  0  from  (2.1).  From  (1.3),  (2.11)  and 
(2.12),  for  Id  <  1  , 


where  I*  now  simplies  to; 


[C  +  C'l 

1 

J 

'  /■/ 

ic'i=i  C 

[c'  -  cJ 

kc(C',0l' 

For  1^1  >  1  ,  (2.17)  simplifies  to: 


(3.2) 


where 


zt  =  qi  Z(  +  92 


9i  (Ci  i)  =  C  h 


(3.3) 

(3.4) 

-  ,1  s  (3.5) 

IT  Zf  (j,() 

Note  that  even  though  has  one  analytic  expression  (3.2),  it  defines  two  different  analytic 
functions  inside  (used  in  (3.1))  and  outside  the  unit  circle  (used  in  91  in  (3.4)). 

^From  the  alternate  expression  (2.16)  for  ,  it  follows  that  for  |(^|  >  1 


92  = - 


2C 


9i(C,0  =  -c  - 


2^‘  (h 


k  (^0 


(3.6) 


Further,  noting  that  ^  and  hence  ^  are  analytic  inside  the  unit  circle  and  t)  and 
•2c(C*5  0  are  complex  conjugates  on  |C|  =  1  ,  it  follows  from  (3.5)  and  Cauchy’s  integral 
formula  that 

or  r  zir>  1 

(3.7) 


a  =  -K  I  <  1 

^  C'-l/C  ^c(C',<) 


^From  (3.2),  (3.4)  and  (3.7)  it  is  clear  that  small  changes  of  ^  on  |(^|  =  1  will 
induce  only  small  changes  in  91  and  92  for  |C|  >  1  and  so  determination  of  91  and 
92  from  ^  on  the  unit  circle  is  well  posed,  when  jz^|  on  |C|  =  1  is  bounded  away  from  zero. 
Further,  note  that  if  we  know  ^(C,^)  a  circle  of  radius  R  where  R  >  1  (Note  R  could 
depend  on  t  )  but  smaller  than  the  distance  of  the  nearest  singularity  of  z^(^,t)  other 
than  that  at  the  origin,  than  for  any  |C|  <  72,  including  |C1  =  1, 


KC,()  =  -  §  I"  <  +  5^  /  77^  WC'.O  +  ;  c'l 

TT  2m  y|^/|  =  R  C  ~  s  ^ 


(3.8) 


for  the  channel  geometry  and  for  the  radial  geometry 

a(t) 


(3.9) 
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where  a{t)  is  determined  from: 


“(<)  =  5^  /  <*<’  <C,  t)  (3.10) 

The  determination  of  z  for  |C|  <  R  from  given  z  on  |C|  =  i2  is  again  well  posed 
from  the  above  formulae. 

Further,  note  from  (3.2)  that  for  |C|  <  1  ,  Re  li  defines  a  harmonic  function  so 
that  on  1(^1  =  1  , 

RehU,t)  = - - - ^  -  (3.11) 

kc(C,0l' 

Thus  from  the  maximum  principle  for  harmonic  function  for  |Cl  <  1  : 

2  1  2  1 

max  —  - o  >  —Re  I4  >  min  —  - 5-  >  0  (3.12) 

ICI  =  1  TT  ^  1  TT 

Further,  a  little  manipulation  of  (3.2)  gives: 


r  1  i  ^  rci  +  C 

■•'o’  ^  2^n  c  [c  -  Cl. 


(3.13) 


It  is  clear  that  Re  l4,{-^,i)  is  a  harmonic  function  in  lCi|  <  1  talcing  on  the  boundary 
value: 

fe/,(i-,t)  =  ^  (3,14) 


on  ICil  =  1  •  By  putting  =  i  ,  from  maximum  principle,  it  follows  that  any 

ICl  >  1 , 

12  .12 

max  - 7  —>  Re  I4  >  min  - 5-  —  >  0  (3.15) 

ICl  =  1  TT  ICI  =  1  TT 


The  properties  (3.12)  and  (3.15)  will  be  crucial  as  we  shall  now  see. 

The  ensuing  discussion  will  focus  on  solutions  to  (3.3)  with  initial  data  obtained  by 
analytically  continuing  physically  relevant  initial  data  z(C,0),  subject  to  the  decomposition 
(1.1)  (or  1.5  for  the  radial  case))  and  the  requirement  that  Zf(C,0)  ^  0  in  jt^l  <  1.  There 
is  no  apriori  reason  for  a  specific  solution  of  the  analytically  continued  equation  (3.3)  for 
t  >  0  to  be  identical  to  the  solution  to  the  original  equations  (satisfying  boundary  condition 
(1.10))  with  decomposition  (1.1)  (or  (1.5)  for  the  radial  case).  Therefore,  as  part  of  our 
discussion,  we  will  include  analytical  arguments  that  suggest  that  the  solution  to  (3.3) 
that  we  calculate  in  |^|  >  1,  on  analytic  continuation  back  to  |^|  <  1,  is  consistent  with 
decomposition  (1.1)  (or  1.5  for  the  radial  case)  and  that  with  a  finite  speed  of  any  zero  of 
there  will  be  no  zero  of  in  |C|  <  1,  atleast  for  a  finite  nonzero  time  interval. 
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Equation  (3.3)  is  in  some  sense  a  nonlocal  hyperbolic  p.d.e,  if  indeed  hyperbolicity  can 
be  defined  for  such  equations.  The  coefficients  are  clearly  analytic  for  |C1  >  1,  except  at 
infinity,  where  from  (3.5)  and  (3.6),  qi  is  proportional  to  (  ,  but  52  is  analytic.  In  reality, 
(3.3)  is  a  nonlinear  integro-differential  equation.  The  analyticity  of  the  coefficients  qi  and 
52  >  however,  guarantee  that  there  is  no  spontaneous  generation  of  singularity  in  |C|  >  1  for 
if  it  were  otherwise,  reversal  of  time  would  imply  that  the  singularity  moves  somewhere  else 
rather  than  disappear  (see  Lacey,  1982  for  a  discussion  of  moving  singularities  in  a  differing 
but  equivalent  formulation).  Thus  any  singularity  of  z(C,t)  outside  the  unit  circle  at  a 
point  must  be  the  result  of  singularities  moving  from  somewhere  else.  Earlier,  Richardson 
(1972)  showed  the  invariance  of  singularity  form  using  a  Schwartz  function  approach.  Even 
though,  he  was  concerned  with  a  displacement  of  a  less  viscous  fluid  by  a  more  viscous 
fluid  rather  than  the  other  way  around,  his  results  hold  by  using  time  reversal  argument. 
However,  the  motion  of  singularities  was  not  addressed  in  the  Richardson  (1972)  paper. 

We  note  from  (3.3)  that  the  characteristic  velocity  is  —qi  that  is  always  pointed 
inwards  towards  |(|  =  1  . 

^From  (3.4)  and  (3.15),  it  follows  that 


Re 


Ut) 

Cc(<) 


=  -  Re  UiCc{t),t)  <  0 


(3.16) 


on  a  characteristic  (c{t)  -  Equation  (3.16)  immediately  implies  that  the  radial  component 
of  cheiracteristic  speed  at  any  time  t  at  any  ^  outside  the  unit  circle  is  always  inwards 
towards  |C|  =  1  and  so  information  in  |C|  >  1  flows  inwards  towards  |C|  =  1  at  all 
times  for  which  a  solution  exists.  In  peirticular,  this  implies  that  singularities  of  the  initial 
conditions  in  |C|  >  1  move  towards  the  unit  circle  since  a  singularity  must  move  with 
the  characteristic  speed  —qi  evaluated  at  that  point. 

On  the  other  hand,  for  |^|  <  1  ,  for  given  ,  (3.1)  is  a  hyperbolic  equation  with 
characteristic  speed  — C  I4  ,  which  according  to  (3.12),  moves  outwards  towards  1C|  =  1  • 
Further  the  characteristic  speed  is  zero  at  ^  =  0  so  any  singularity  of  z  at 
C  =  0  at  initial  time  will  not  move  in  time. 

^From  the  exterior  equation  (3.3),  as  C  ^  00  >  using  (1.1),  (3.5)  and  (3.6)  the 
asymptotic  equations  are: 

=  5  (Z;  +  2  (3.17) 

for  the  channel  case  and  for  the  radial  geometry  from  (1.5),  (3.5)  and  (3.6),  we  get: 


zt 


q'(0 

a{t) 


C  ^<  + 
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4 

TT  a(f)C 


(3.18) 


The  general  solution  to  the  initial  value  problem  in  (3.17)  is 


2(c,<)  =  2(  +  0) 


(3.19) 


and  for  (3.18),  it  is: 

Thus  any  singularity  present  at  infinity  at  initial  time  does  not  move  to  finite  (  values. 

Now  consider  the  motion  of  a  zero  Co(^)  of  ■  Here,  by  a  zero  of  ,  we  only  refer 
to  points  where  Z(;  is  zero  but  analytic.  If  z^  is  zero  but  nonanalytic  at  a  point  (say  like 
a  branch  point),  then  it  is  called  a  singular  point  and  will  be  denoted  as  (^,(<)  just  like 
any  other  singularity. 

To  follow  the  motion  of  zeros,  it  is  convenient  to  take  the  derivative  of  (3.3)  with 
respect  to  (  and  divide  the  resulting  expression  by  —z^  to  obtain: 


(h  -  (h 


iii 


(3.21) 


iFrom  the  form  (3.21),  one  notes  that  a  double  zero  of  z^  is  only  possible  when  92 =  0. 
This  can  occur  momentarily;  however,  we  limit  our  discussion  to  a  simple  zero  for  which 


~  «cc(Co(<)>0  (C  -  Co(<)) 


(3.22) 


On  substituting  this  into  (3.21),  we  find 


Co  =  -  gi(Co(t),t) 


g2c(Co(0»0 

2^<c(Co(0>0 


(3.23) 


We  are  unable  to  show  in  the  general  case  that  Me  ]  is  negative,  i.e.  if  the  zero  must 
always  approach  the  physical  domain.  Also,  another  open  question  is  whether  the  velocity 
of  approach  is  always  smaller  than  —gi(Co(i),  t)  ,  which  would  have  been  the  velocity  of 
approach  of  a  singularity  at  that  point.  However,  this  is  the  case  for  the  specific  case  of 
SaiFman  (1959)  exact  zero  surface  tension  solution  (in  the  formulation  of  this  paper).  In 
that  case 

z(C,  t)  =  i  +  d{t)  -  -  In  C  +  -  {I- \)  In  {I- ai  (t)  C")  (3.24) 

TT  TT 

where  1  >  ai(0)  >  0  and  0  <  A  <  1  ,  and  d{t)  and  ai(t)  are  determined  by 


X  d  4"  —  (1  —  A)  In  ai  —  t  Kq 

TT 


(3.25) 
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(3.26) 


d  —  —  In  ai  +  —  (1  —  A  )  /n  (1  —  aj)  =  Ki 

TT  TT 

where  Kq  and  Ki  and  depend  on  the  initial  conditions.  In  this  case,  for  0  <  A  <  A  , 
we  calculate  Co  =  ±  *  on  the  imaginary  axis  and  corresponding  speed  Co  at  all 

V  ai(l— 2A) 

times  is  towards  to  the  physical  domain  |C|  =  1  though  the  magnitude  of  the  speed  is 
smaller  than  the  magnitude  of  — 9i(Co(t),t)  ,  calculated  from  (3.6)  which  is  also  directed 
towards  the  physical  domain  |Cl  =  1  .  For  1  >  A  >  ^  ,  within  the  Saffman  class  of 


solutions,  Co  =  ±- 


,  on  the  real  positive  axis  at  a  distance  further  from  any  of 


’  ”  - ” 

the  poles  at  C  =  Cs  =  speed  of  each  zero  is  directed  towards  unit  circle, 

though  they  are  less  than  |  —  qi  (Co(t))  t)!  in  absolute  value.  Within  this  class  of  solutions, 
as  t  oo,  the  zeros  actually  eventually  settle  down  at  C  =  ±  v/2A-i  ’  finite  distances 
from  the  physical  domain  |C|  =  1  .  This  is  unlike  other  cases  given  by  Lacey  (1982), 
Howison  (1985,1986ab),  Shraiman  Sz  Bensimon  (1985)  where  the  zeros  actually  impinge 
the  physical  domain  in  finite  time. 

Now  consider  the  initial  value  problem 


^(C^O)  = - ln(  +  i  +  /„(0)  C” 

TT 


(3.27) 


for  the  channel  geometry  or 


<c,o)  =  ^  +  x; 

^  n=0 


(3.28) 


for  the  radial  geometry.  We  assume  that  2^  is  nonzero  for  |Cl  <  Ro  and  that  the 
convergence  of  the  infinite  series  occurs  within  a  circle  of  radius  Rq  (  >  1).  When  this 
series  is  not  convergent,  we  can  take  its  unique  einalytic  continuation  across  the  unit  circle 
to  be  the  initial  condition.  As  long  as  integrals  (3.2)  and  (3.7)  exist,  each  of  91  and 
$2  are  analytic  functions  outside  the  unit  circle  except  at  C  =  00  where  qi  is 

proportional  to  ^  (92  is  analytic).  Thus,  any  singularity  for  2  outside  the  unit 

circle  has  to  be  there  at  initial  times  and  owing  to  the  hyperbolic  nature  of  the  equation, 
any  such  singularity  at  (s{t)  must  move  with  speed  ,  which  owing  to  the 

property  (3.16)  is  always  towards  the  unit  circle.  This  is  the  basic  reason  for  ill-posedness 
for  the  initial  value  problem  restricted  to  [Cl  <  1  since  two  slightly  different  data  in 
ICI  ^  1  can  correspond  to  different  initial  distribution  of  singularities:  in  one  case  the 
singularity  may  be  far  off  while  in  the  other  case,  the  singularity  can  be  made  as  close 
as  we  want  to  the  physical  domain.  After  a  short  time,  in  the  latter  case,  depending  on 
the  singularity  nature,  it  can  impinge  the  physical  domain  and  the  difference  between  the 
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interfacial  shapes  in  the  two  situations  in  any  norm  that  measures  the  interfacial  slope  will 
be  0(1)  .  The  singularity  location  and  strength  can  be  arranged  so  that  this  happens  for 
arbitrarily  short  times.  However,  when  one  chooses  an  norm  that  distinguishes  between 
initial  conditions  based  on  their  values  in  a  strip  adjoining  the  physical  domain  then  up  to 
certain  times,  the  evolution  problem  can  be  expected  to  be  well  posed. 

However,  the  problem  is  well  posed  when  |C|  >  1  is  our  domain  as  we  now  argue. 
We  assume  that  there  exists  T  >  0  so  that  the  speed  of  a  singularity  (given  —qi  )  or 
a  zero  (given  by  3.22))  has  uniform  upper  bound  M  for  the  inward  radial  component  for 
0  <  t  <  T  .  This  appears  to  be  reasonable  assumption  since  all  the  solutions  known 
explicitly  satisfy  this  criteria.  Then  if  we  choose  Ti  <  T  so  that  Rq  —  MT\  >  1  ,  then  it 
is  clear  that  for  te[0,  Ti],  Ro—Mt  >  1  so  that  z  is  analytic  for  1  <  |C|  <  Ro  —  M  t . 
For  such  a  time,  each  of  qi  and  q^  is  determined  by  first  determining  z  and  therefore 
z^  on  1^1  =  1  using  (3.8)  (or  (3.9)  and  (3.10)  for  the  radial  case)  and  then  using  (3.4)  and 
(3.7).  Further,  the  value  of  z  on  a  circle  of  radius  Rq  ~  M  t  will  completely  determine 
z  in  any  other  circle  with  radius  between  1  and  Rq  —  M  t  through  (3.8)  (or  (3.9)  and 
(3.10)  for  the  radial  flow  ).  This  determination  is  well  posed  as  argued  earlier.  Further,  it 
is  clear  that  the  values  of  z  at  time  t  determine  z  at  any  points  on  the  closed  curve 
5  at  time  i  shown  in  Fig.  5.  The  curve  S  is  chosen  such  that  it  determines  z  on  the 
circle  of  radius  Rq  —  M  (t  +  6t)  ai  time  t  +  6  t  by  using  the  method  of  characteristic 
on  (3.3).  This  determination  is  well  posed  since  a  characteristic  method  is  equivalent  to 
an  ordinary  differential  equation,  for  which  continuous  dependence  on  initial  values  is  well 
established. 

The  arguments  in  the  last  paragraph  suggests  the  well  posedness  of  the  formulation 
of  the  exterior  problem.  However,  one  must  ascertain  that  this  exterior  solution  to  the 
analytically  continued  equation  is  actually  in  the  form  (1.1)  or  (1.5)  for  Ti  >  t  >  0  and 
z^  has  no  zero  in  1^1  <  1  as  must  be  the  case  for  a  physically  acceptable  solution. 
This  follows  from  noting  that  (3.1)  is  actually  the  analytic  continuation  of  (3.3)  back  into 
1^1  <  1  across  |C|  =  1.  Further,  at  initial  time,  by  the  very  construction  of  z((,0),  the  data 
satisfied  in  the  interior  problem  is  consistent  with  (3.27)  (or  (3.28)  for  the  radial  flow  case) 
and  as  we  mentioned  before,  the  characteristics  in  the  interior  problem  is  always  pointing 
outwards  with  zero  characteristic  speed  at  the  origin.  This  means  the  solution  of  (3.1) 
cannot  develop  singularity  unless  I4  is  not  analytic  or  fails  to  exist  for  |C|  <  1.  However 
from  the  exterior  problem  ^  is  analytic  upto  time  Ti  on  jCl  =  1  and  in  particular 
must  be  integrable  implying  I4  exists  and  defines  an  analytic  function  in  |C|  <  1  . 

Thus,  in  the  interior  evolution  the  only  singularity  is  at  the  origin  in  the  form  given  in 
(1.1)  (or  (1.5)  for  the  radial  geometry).  Further,  no  zero  of  z^  can  form  in  |C|  <  1  in  the 
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time  of  interest  since  at  initial  time  there  is  no  such  zero  and  the  only  way  for  a  zero  to 
form  in  this  domain  is  for  2!ri  ijc'l  =  1  ^  C'  to  change  values  discontinuously  in 

time.  This  can  happen  only  if  a  zero  of  moves  through  |Cl  =  1-  Clearly,  none  of  the 
above  is  the  the  case  since  from  the  exterior  problem  and  its  inverse  are  analytic  on 
the  unit  circle  for  <  e  [0,  Ti]  .  Thus  the  solution  one  calculates  in  the  exterior  is  actually 
of  the  form 

2  °° 

z(C,i)  =  —  (n  C  +  i  +  V  /„(<)  C"  (3.29) 

for  the  channel  geometry  and  of  the  form 

.=(<.<)  =  ^  +  f;  M«)C"  (3.30) 

for  the  radial  geometry  for  |^|  <  1  with  Z(^  ^  0  in  that  domciin.  This  guarantees  that 
the  initial  value  problem  solved  in  |(^|  >  1  ,  when  analytically  continued  to  |^|  <  1  ,  is 
in  accordance  to  the  requirements  on  in  |C|  <  1. 

As  a  concrete  example  of  how  one  might  solve  an  initial  vzilue  problem  with  branch 
point  singularities,  we  start  with  initial  value 

=  E  Trrr; «  -  0(o))‘'^'  £j(c,o)  +  g(c,o)  (3.31) 

where  |Cj  |  >  1  and  l3j  is  generally  a  noninteger  constant  and  £?j(Cj(0),0)  ^  0.  In  the 
case,  when  /3j  =  1  ,  we  replace  the  above  expression  by  the  limit  of  — v  1  which  clearly 
gives  a  logarithmic  expression.  Further,  suppose  £?j((^,0)  is  entire  and  G(C,0)  is  an 
analytic  function  everywhere  in  the  finite  C  plane  except  at  C  =  0  ?  where  the  singularity 
is  consistent  with  the  singularity  of  5:  as  given  by  (1.1)  (or  (1.5)).  On  substituting 

j=i  ^ 

into  (3.1),  where  satisfies 

0  =  -«i(Cy(0.0  (3.33) 


and  each  Ej((^,t)  is  chosen  to  satisfy 


p  ^  p  _  n  (?'!  F  “  9i(Cj(<),<)] 

-  qiEj^  -  - 


(3.34) 


then,  it  becomes  clear  that 


Gt  —  qi  G(^  +  q2 


(3.35) 
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On  examination  of  (3.34)  and  (3.35),  it  is  clear  that  for  times  when  \Cj{t)\  >  1  ,  and 
^  0  in  1^1  <  1  ,  each  of  Ej  and  G  will  remain  analytic  in  |^|  >  1  since  51  and 
52  are  remain  analytic  outside  the  unit  circle. 

We  now  outline  how  a  rigorous  proof  of  the  analyticity  would  follow.  Consider  G(C,  i)- 
If  dC  G((^,t)  ^  0  for  some  finite  contour  C  completely  outside  jC|  >  1,  the  preimage 
Co  (which  must  be  a  closed  contour)  of  the  contour  C  under  the  flow  (i-e.  each 

point  of  the  contour  moved  according  to  the  equation  (  =  —gi(Ci  i))  time  t  =  0  must 
also  be  a  flnite  contour  entirely  in  |<^|  >  1  by  noting  the  property  that  Re  [gi/C]  >  0  and 
finite  for  the  time  that  a  zero  or  a  singularity  of  z^  has  not  impinged  |C1  =  1-  However, 
d(  G((,  0)  =  0  from  the  assumed  analyticity  initially.  On  the  otherhand  it  is  clear 
from  (3.35)  and  the  analyticity  of  92  in  ICl  >  1  that  G(C,  i)  =  0  on  a  contour  moving 
with  the  flow  — gi(C>f)-  This  leads  to  a  contradiction,  implying  that  G(C,t)dC  =  0 
for  arbitrary  closed  contour  C  in  |^|  >  1,  which  by  Moerara’s  theorem  implies  G{^,t)  is 
analytic  in  |Ci  >  1.  The  proof  for  each  Ej  would  be  along  similar  lines. 

Thus  the  form  of  the  solution  (3.32)  implies  invariance  of  the  singularity  form  at  later 
times.  Note  however,  that  inside  the  unit  circle  at  ^  ^  »  ^^ch  of  G  and  Ej  may 

have  have  singularities  for  t  >  0  ,  as  the  analytic  continuation  of  gi  and  §2  has 
singularities  inside  (from  (3.5)  and  (3.6)).  However,  since  z(C,t)  as  argued  earlier  cannot 
have  singularity  inside  except  at  (^  =  0  ,  the  singuleirities  of  Ej  ’s  and  G  must  cancel 
out  at  C  =  ^  on  the  specific  Riemann  sheet  corresponding  to  the  physical  domain. 

Given  Ej  and  G  at  any  instant  of  time  on  a  circle  of  radius  R  greater  than  unity, 
but  smaller  than  the  nearest  zero  or  singularity  of  z^  ,  (3.32)  determines  z  on  that  circle 
which  then  determines  z  and  hence  z^  on  (C|  =  1  through  (3.8)  (or  (3.9)  and  (3.10)  for 
the  radial  geometry).  z<;  on  |C|  =  1  determines  gi  and  92  from  (3.2),  (3.4)  and  (3.7). 
Equations  (3.33)-(3.35)  can  then  be  used  to  advance  Ej  and  G  along  a  characteristic. 

We  now  discuss  the  fate  of  singularities.  First,  we  assert  that  a  singularity  cannot  stay 
away  at  a  finite  distance  from  |^|  =  1  since  this  would  otherwise  imply  existence  of  a 
limit  point  or  a  limit  cycle  where  Re  [^i  /^]  =  Re  I4  =  0.  From  maximum  principle  and 
(3.15),  it  follows  that  this  is  impossible.  Thus  every  singularity  must  indefinitely  approach 
1^1  =  1  .  We  know  from  Howison  (1985,1986ab)  within  the  context  of  exact  solutions 
that  pole  singularities  of  z^  never  actually  impinge  |(|  =  1  in  finite  time  though  it 
approaches  it  exponentially  in  time.  We  now  discuss  in  a  general  context  if  a  singularity 
Ca(t)  >  once  very  close  to  the  unit  circle,  can  actually  impinge  the  unit  circle  in  finite  time. 
Consider  singularity  in  the  form; 


~  B„(()  (c  -  c.wr'’ 


(3.36) 
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for  C  substitution  into  (3.3),  we  can  easily  establish 

Bo{t)  =  .Bo(O)  -fo  (3,37) 

Note  that  for  a  simple  pole,  i.e.  for  /9  =  1  ,  the  strength  5o(t)  is  a  constant.  This  is 
known  for  the  previous  SHSB  exact  solutions. 

Assume  that  at  time  to  >  the  singularity  is  very  close  to  the  unit  circle.  We  will  now 
find  the  asymptotic  equation  of  motion  for  the  singularity  beyond  that  time.  Let 


Cs  =  Re' 


(3.38) 


Since  singularity  moves  wich  the  characterestics,  C*  =  —qiiCa,t)  ,  which  from  (3.2)  and 
(3.4)  imply  that 

^  +  iu,  =  (3.39) 

We  now  consider  the  asymptotic  expression  for  I4  as  R  1"^  .  ^From  (3.2),  we 
immediately  get 

1  \  R  A-  1  1 

-/,(«<),.)  =  ^  I  ..  (3  «) 

Let’s  assume  ^  ±7r  or  otherwise  the  range  of  integration  on  u  will  be  changed 
from  0  to  27r  and  the  same  argument  as  given  below  can  be  repeated.  We  break  up  the 
integral  in  (3.40) 

/TT  ev,—f.  rir  /•!'•+« 

=  /  +  /  +  /  (3.41) 

-TT  J —n  Jv,  +  t  J  v,—t 


where 


€  =  (iJ-1)’/^ 


(3.42) 


It  is  then  clear  that  R—1  «  e  <<1  and  the  integrand  in  each  of  the  1st  two  integrals 
can  be  replaced  by 

-icod(u  -  i/,)  -  hR  -  1)  cosec^hu  -  u,)]  2  +  o{R-l)  (3.43) 

On  integration  by  parts  of  the  second  term  in  (3.43),  one  finds  that  the  asymptotic  behavior 
for  ^  >  I  of  the  sum  of  the  first  two  integrals  in  (3.41)  is  given  by 

/"•“*/*  i  \  1  1 

■•■  / - ^  /  ducot-(u-v,)  —--.2  +  o(l) 

-JT  Jv,+f  ^  J —‘K  2  |z^(e*  ,t)| 
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1 


+  o(J?-l)  (3.44) 


(j^-1)  r 


dv  Cot-{u  -  Ua) 


d 

dv  lz^(e«-,Ol" 


The  first  term  above  is  completely  imaginary  and  will  not  be  relevant  to  our  discussion  of 
the  asymptotic  equation  for  R  found  by  taJcing  the  real  part  in  (3.39).  For  <  |  , 
one  finds  from  integration  of  (3.43)  the  real  part  of  of  the  first  two  integrals  in  (3.41) 
is  o((i?  —  1)^^)  .  Now  consider  the  contribution  from  the  third  integral.  Introduce  new 
variable  s  =  ,  we  find  that  to  the  leading  order 


+  e 

c 


2  (fl  -  1)^'’  /•’*f  .  (-!-»)  ,,  ,  2 

:r2  |B„(()|^  (l+»^)  + 


2  cos(s[i2  —  1]) . 


(3.«) 


For  3  <  1 


to  the  leading  order  the  above  expression  simplifies  to 


\Bo{t)\^  Jo 


^-1 


(3.46) 


On  the  other  hand  for  ^  >  i  ,  it  is  clear  from  (3.45),  that  the  contribution  is 

0{{R—\)^)  which  is  smaller  than  0(R  —  1)  contribution  to  Re  in  (3.44).  Combining 
the  above  information  with  the  real  part  of  (3.39)  we  obtain 


R - Mi{t){R-lf^ 


(3.47) 


for  J?  ^  I"*"  for  /3  <  I  where 

4  1 

=  (1+*)"  (3.48) 

|Bo(«)|  Jo 

and  for  )9  >  ^  , 

it - M2{t)  {R  -  1)  (3.49) 

where 

Mm  =  1  d.  cotl  (.  -  ..)  (3.50) 

Note  that  from  (3.15)  it  follows  that  M2{t)  as  determined  above  must  be  positive. 
Further  from  (3.48),  for  any  v(t)  ,  Mi{t)  cannot  go  to  zero  when  Bo{t)  stays  bounded. 
Further  from  (3.37),  and  the  relation  of  with  li  through  (3.4),  it  follows  that  Bo{t)  will 
stay  bounded  when  ^  >  0  .  For  ^  <  0  ,  it  is  possible  to  have  So(t)  unbounded  as 
— ♦  1.  y,From  (3.50),  it  is  clear  that  that  will  stay  bounded  for  M2{t)  when  /3  >  |  . 
This  is  true  irrespective  of  the  dependence  of  Ug  on  t  . 
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On  solving  these  ordinary  differential  equation,  we  find 


R-1  =  [-{1-20)  f  Mi{t')dt'  +  (il(0)  -  1)'-"^] 

Jto 

for  yd  <  I  and  for  yd  >  |  ,  we  find 


(3.51) 


R  -  I  =  {Ro-l)e  ^‘0  (3.52) 

It  is  clear  that  for  0  <  ^  ,  ii  —  1  will  be  zero  unless  (C«(OiO  — »  oo  eis 

|C«|  —*  1  such  that  Bo(t)  as  determined  from  (3.37)  -+  cx>  in  such  a  way  that  the 
integral  of  Mi  in  (3.51)  stops  growing  indefinitely  as  t  — »  oo  .  We  cannot  rule  out 
such  an  exceptional  case  when  0  <  0  ,  since  the  behavior  of  /4((^,(t),t)  given  by  (3.46) 
implies  that  9io,^  (Ca(0?0  determined  from  (3.4)  will  tend  to  infinity  as  |(a|  —*  1  for 
yd  <  0.  However,  more  careful  arguments  not  considered  here,  may  rule  that  out  as  well 
though  we  are  unsure  about  it.  Thus  aside  from  the  exceptional  case  pointed  out  over 
here,  there  will  be  finite  time  singularity  {ox  0  <  |  . 

^From  (3.52),  for  /3  >  i  the  singularity  never  hits  the  the  physical  domain  in  finite 
time  as  M2(<)  in  (3.50)  is  always  finite.  At  the  cross  over  value  of  yd  =  ^  ,  the  three 
integrzd  terms  contribute  equally  and  one  finds  a  (i?  —  1)  contribution  for  which  (3.40) 
holds  (though  with  a  different  M  ),  and  therefore  in  this  case  as  well,  the  singularity  will 
not  reach  the  physiczd  domain. 

In  the  case  of  ^  <  0  ,  if 

HiM  ~  +  Bo(t)  (C-C. (<))■'  (3.53) 

as  C  C«(0  nonzero  EQ{i)  ,  then  it  is  easy  to  see  that  Re  l4((t(t),t)  =  0(1)  . 

From  the  real  part  of  (3.39),  it  follows  that  R  —  1  in  this  ceise  will  be  zero  zdter  a  finite 
time.  Thus  for  a  singularity  of  the  form  (3.53)  with  0  <  0  ,  it  will  move  to  the  boundary 
of  the  physical  domain  in  finite  time. 

We  now  sketch  a  numerical  algorithm  to  solve  for  z{(^,t)  outside  the  unit  circle. 
We  will  limit  our  description  to  the  channel  geometry  though  the  procedure  with  minor 
chtinges  is  readily  applicable  to  the  radial  flow.  We  will  also  only  describe  the  algorithm 
to  calculate  z{(,t)  on  a  circle  of  radius  R{t)  ,  where  R{t)  shrinks  with  t  in  some 
prescibed  way  and  initially  does  not  contain  any  singuleirity  of  Z(^  or  a  zero.  We  will  find 
that  we  can  stably  calculate  z{^,i)  on  the  unit  circle,  i.e.  the  physical  interface  as  a  by 
product. 
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Assume  N  uniformly  spaced  out  points  in  a  circle  of  radius  R{t)  at  any  time  at 
which  is  known.  We  now  describe  how  to  calculate  z{C,t  +  S  t)  on  a  circle  of 

radius  R{t  +  6  t)  ,  where  R{t)  is  changed  in  some  prescibed  way. 

^PVom  a  discrete  fast  Fourier  transform,  we  can  calculate  each  cj  in  the  relation: 

2  ^ 

z  +  -  In  (  =  ^  Cj  W  e‘>  "  (3.54) 

^  —N  I 

where  ^  =  Re""  at  points  on  the  circle  of  radius  R{t)  .  We  know  from  previous  discussion 
for  the  appropriate  initial  value  problem  for  the  continuum  case  that  cj  =  0  for  j  <  0  . 
This  should  be  approximately  true  for  the  discrete  version  for  large  JV  .  This  provides  a 
check  on  numerical  accuracy.  On  finding  the  cj  ,  we  calculate  z  and  for  =  e'"  on 
N  uniformly  spaced  out  points  using 


2  * 

2  +  -  fn  C  =  5];  C,  e«  >  "  (3.55) 

^  >=o 

2  1  ^ 

+  -7  =  12  e‘  ^  (3.56) 

^  i=o 

Note  that  the  Cj  for  negative  j  is  not  used.  We  now  can  calculate  92  at  any  point 
outside  the  unit  circle  using  (3.5)  as  2  =  i  for  the  channel  geometry  and  therefore 


92 


C 


-I  c  +  E7=o  j  cj 


(3.57) 


Further,  on  finding  2^  and  hence  ^  at  iV  uniformly  spaced  out  points  on  the  unit 
circle,  we  \Jculate  a  finite  set  of  coefficients  dj  and  6j  (  iV  of  them  altogether)  in  the 
Fourier  series  representation  of 


2 

7r|2^(e*‘',<)|^ 


00 

do  +  ^  dj  cos  (j  1/  +  (5j) 
>=i 


^From  (3.2)  and  (3.4),  it  follows  that 


9i  =  -C 


do  +  d,  C"-' 

j=i 


(3.58) 


(3.59) 


After  knowing  dj  in  (3.58),  a  finite  truncation  of  (3.58)  can  be  used  to  calculate  q\  at 
any  C  point  outside  the  unit  circle. 
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Now,  we  put  (3.3)  in  a  characteristic  form  and  consider  the  problem  of  determining 
z  on  a  circle  of  radius  R{t  +  6  t)  chosen  such  that 

—In  R{t  +  6t)  +  R(t)  =  max  Re  ^  6t  (3.60) 

Id  =  R{t)  L  C  J 

Now  the  val  ue  of  z  at  a  point  P  at  time  t  S  t  in  Fig.  5  is  determined  from  value 
of  2  ,  gi  and  92  at  point  P\  on  the  curve  S  at  time  t  and  qi  ,  52  at  P  at  time 
t  +  St  .  Applying  a  second  order  difference  scheme  for  the  characterestic,  we  obtain: 

C(P)  -  C(Pi)  =  -5  l9i«(A).i)  +  9i(C(P),  t  +  «))  «  (3.61) 

2(C(P),<  +  «)  =  2«(Pi),()  +  I  b(C(Pi),()  +  92(C(P).*  +  «)1  St  (3.62) 

In  the  above,  C(Fi)  Eind  z(P,t  +  St)  for  N  uniformly  spaced  out  points  P  on  the 
circle  of  radius  R{t  St)  have  to  be  determined  all  at  once  in  an  iterative  fashion  since 
each  of  qi  and  92  at  time  t  +  St  are  determinable  in  terms  of  2  at  t  +  St  on 
N  uniformly  spaced  out  points  on  circle  of  radius  R(t  +  St)  (as  described  in  the  last 
peiragraph). 

The  procedure  can  be  repeated  until  R{t  +  St)  <  I  ov  Z(^  =  0  on  the  unit 
circle  for  which  qi  will  fail  to  exist.  We  note  that  to  follow  the  motion  any  singularity, 
it  is  not  necessary  to  calculate  2  near  that  singularity.  Once  qi  is  calculated  as 
in  the  above  procedure,  it  is  known  for  every  ^  outside  the  unit  circle  and  we  can 
then  cjilculate  each  singularity  location  by  numerically  integrating  the  ordinary  differentiid 
equation  =  — 5i(Ca(0j0  •  The  strength  of  the  singularity  Bo{t)  given  by  (3.37)  can 
similarly  be  calculated  at  any  time  through  integration. 

In  the  following  sections,  we  consider  how  small  nonzero  surface  tension  modifies  the 
dyneunics  of  z{C,t).  Our  discussion  of  surface  tension  effects  will  be  limited  to  initial 
conditions  that  are  independent  of  B  .  Further,  we  will  assume  that  the  singularity  of 
z^  at  initial  time  is  either  a  simple  pole  or  a  branch  point  of  a  power  type  like  any  of 
the  ill  the  expression  (3.31)  with  0  <  <  1.  We  will  eilso  assume  that  the  zeros  and 

singularities  of  zq(^  are  at  0(1)  distance  from  each  other.  Further,  our  results  in  sections  4 
to  7  are  restricted  to  cases  when  the  distance  of  each  singularity  of  Z(^  from  the  boundtiry 
of  the  physictil  domain,  K]  =  1  ,  is  0(1).  However,  since  the  physical  interface  is  |Cl  =  1, 
the  dynamics  described  in  these  sections  are  not  directly  relevant  at  that  instant  of  time  to 
the  large  distortions  of  interfacial  shape  caused  by  one  or  more  singularities  coming  very 
close  to  Id  =  1.  For  large  times,  all  singularities  must  come  very  close  to  (d  =  1  and  the 
study  of  the  this  latter  stage  is  taken  up  in  section  8,  though  in  a  limited  way. 
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4.  Perturbation  expansion  in  powers  of  B  . 

The  discussion  in  the  last  section  suggest  that  the  zero  surface  tension  problem  solved 
in  the  unphysical  plane  is  well  posed.  Thus,  it  is  natural  to  study  the  effect  of  small  non 
zero  surface  tension  in  this  unphysical  plane.  We  start  by  assuming  that  ein  asymptotic 
expansion  for  some  part  of  the  complex  ^  plane  at  least  in  the  vicinity  of  the  physical 
domain  is  of  the  form 

z  zq  B  zi  +  ...  (4.1) 

In  this  section,  we  will  examine  where  such  a  perturbation  expansion  is  inconsistent.  In 
latter  sections,  we  then  proceed  ^’o  show  that  under  certain  assumptions,  a  so  called  ’inner’ 
region  exists  (in  what  is  now  the  standard  language  in  the  theory  of  matched  asymptotic 
expansion)  such  that  an  appropriate  inner-outer  matching  is  possible. 

On  substituting  (4.1)  into  (2.17),  it  is  clear  that  zq  denotes  the  z  determined 
for  B  =  0  in  the  last  section.  Going  back  to  (2.17),  it  is  clear  that  zi  satisfies  the 
integro-differential  equation: 


^1*  “  9io  —  52i  +  9li  2:o<;-|-93o  + 


^4o 


950 


3 

2  z^^^n 


(4.2) 


where  subscript  0  on  each  of  ,  53  ,  54  ,  and  qj  denote  the  evaluation  of  this 
quantities  as  given  by  (2.18),  (2.20)-(2.23)  using  =  0  and  z  =  20  as  determined 
in  the  last  section.  The  subscript  1  on  qi  and  52  denotes  the  derivatives  of  each  of 
qi  and  92  with  respect  to  5  at  5  =  0  as  determined  from  (2.18)  and  (2.19)  on 
substituting  (4.1)  for  2  .  We  can  solve  (4.2)  in  the  unphysical  plane  by  a  variation  of  the 
numerical  method  outlined  in  the  last  section  for  zero  surface  tension  since  the  operator 
on  the  left  hand  side  is  similar  to  that  of  the  equation  for  zq  (in  3.3).  It  is  clear  from 
(4.2)  that  at  points  where  20^  is  either  singular  or  has  a  zero,  there  is  singular  forcing  on 
the  right  hand  side  of  (4.2).  In  the  case  when  20^  BMO  (C  -  C.(*))'®  as  (  ^  C.(()  , 
with  (3  <  ^  where  Bo{t)  =  fjj((^a(<),t)  for  determined  by  (3.34)  with  the 

identification  =  C«  >  we  find  from  (4.2) 


^4.3) 


where  Bi{t)  is  determined  by  the  first  order  differential  equation 

=  j97,(WO.<)Bo(<r''"/3(/S  +  2)(/3  +  4)  (4.4) 

with  initial  condition: 

5i(0)  =  0  (4.5) 
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Note  that  ^5  =  1  ,  i.e.  the  singularity  is  a  simple  pole  for  2(;  ,  Bo(t)  is  a  constant  as 
noted  before  in  the  context  of  SHSB  solutions.  For  |  <  ^  <  2  ,  the  most  singular 
term  on  the  right  hand  side  of  (4.2)  comes  from  gij  zq^  rather  than  the  curvature  terms. 
However,  the  uniformity  of  (4.1)  arising  because  of  such  singular  terms  is  simply  because 
the  singularity  of  the  outer  expansion  zq  B  zi  travels  with  speed  qi  rather  than 

gio  .  Thus  in  approximating  the  location  of  the  singularity  by  ,  we  commit  an  error 
in  the  location  of  the  singularity  of  order  B  ,  which  shows  up  in  the  regular  expansion 

(4.1)  in  term  zi  as  a  more  singular  term  ((^  —  ^4)  than  that  of  •  For  /9  in 

the  interval  (4/3,  2)  ,  this  happens  to  be  more  singular  than  the  curvature  term  which 
also  induces  in  zi^  a  more  singular  term  than  Zq^  .  This  nonuniformity  can  be  remedied 
in  the  two  term  outer  perturbation  zq  +  B  zj  by  determining  the  singularity  position 
(s(t)  of  Zq  using 

6  =  -«u«.(<).0  -  (4.5) 

However,  even  with  this  remedy,  there  are  still  singular  terms  on  the  right  hand  side  of 

(4.2)  arising  from  the  curvature,  the  most  singular  of  which  induces  a  singularity  of  the 

form  (4.3)  which  is  more  singular  than  zo<  when  /3  <  2  .  For  0  >  2  ,  there  is 
no  nonuniformity  of  (4.1)  owing  to  curvature,  as  the  induced  singularity  in  is  less 

singular  than  zq^^  in  that  case. 

The  conclusion  is  that  if  the  initial  value  problem  has  a  singularity  for  z^  ~  .Bo(O)  ((  —  Cg(0)) 
I  (a  (0)1  >  1  ,  and  if  that  initial  condition  is  independent  of  B  ,  then  the  singularity  form 
is  preserved  by  zq^  though  its  strength  BQ{t)  generally  changes  (simple  pole  is  an  excep¬ 
tion)  in  time;  the  higher  order  perturbation  term  has  a  worse  singularity  at  approximately 
the  same  point  (to  order  B  )  for  0  <  2  .  Thus,  one  can  expect  an  inner  region  near  any 
such  Ca(0  )  where  the  asymptotic  expansion  (4.1)  becomes  invalid.  This  inner  region  will 
be  treated  in  the  next  section. 

Going  back  to  (4.2),  we  find  that  the  right  hand  side  forcing  terms  is  singular  when 
zq^  is  analytic  but  zero  (generically  simple).  Thus,  if  we  consider  an  initial  value  prob¬ 
lem  for  z(C,t)  with  initial  conditions  that  are  independent  of  surface  tension  B  such 
that  Z((C,0)  =  0  at  some  point,  then  we  know  from  section  3  that  B  =  0  solu¬ 
tion  zo^  preserves  this  zero  (o(t)  that  move  in  time  according  to  (3.23).  At  such  a 
point  (o(0  5  one  finds  a  singular  forcing  on  the  right  hand  side  of  (4.2).  Near  Co(t)  , 
^0<(C»*)  ~  (Co(05  0  (C  ~  Co(0)  and  this  induces  a  singular  term  in  zj  at  Co(t)  so 

that  the  leading  order  asymptotic  behavior  is 

^ic(C,*)  ~  M‘)  «  -  (4.6) 

as  C  Co(0  obtained  by  substituting  (4.6)  into  (4.2)  and  equating  the  most  singuleir 


term.  On  using  (3.23),  we  find  that 


3  [goa(Co(0.<)]  ^  g7o(Co(0>0 
2  q2c(Co{t),t) 


(4.7) 


However  (4.6)  does  not  satisfy  the  initial  condition  zi(C5  0)  =  0  ,  as  it  must  for  B  in¬ 
dependent  initial  condition  ^(C>0)  •  Therefore,  to  a  particular  solution  that  behaves  like 
(4.6),  one  must  add  in  a  solution  to  the  homogeneous  equation  on  left  hand  side  of  (4.2) 
so  that  initial  condition  is  satisfied.  The  homogeneous  part  of  the  solution  must  have 
a  singularity  at  ^  =  Co(0)  at  t  =  0  so  as  to  cancel  the  singularity  in  (4.7)  at 
t  =  0  .  However,  if  there  is  a  singularity  present  in  the  homogenous  solution  at  t  =  0  , 
it  must  be  move  at  the  singularity  speed  — (from  left  heind  side  of  (4.2))  which  is 
different  from  the  speed  of  a  zero  given  by  (3.23).  Combining  all  the  above  information, 
we  conclude  that  for  t  >  0  ,  there  must  be  a  singular  point  Csi  of  that  initially 
coincides  with  Co(0)  ,  i.e.  Cai(O)  =  Co(0)  ,  whose  location  at  later  times  is  determined 
by  Cai  =  -9io  (Cai(^)iO  such  that  as  C  C*i(f)  , 


~  M‘)  «  -  (4.8) 

where 

Ai(t)  =  Ai(0)  e"^  fo  (4  9) 

and  44i(0)  -t-  i4o(0)  =  0.  Thus  one  observes  the  birth  of  a  new  singular  point  Csi(f)  of 
the  outer  asymptotic  expansion  (4.1)  that  is  neither  a  singular  point  of  zp  nor  a  zero  of 
zo(  ,  except  at  initial  time.  Thus,  aside  from  each  singular  point  of  zq^  discussed  earlier, 
an  inner  region  is  required  near  each  zero  Co(t)  zero  of  zqc  (which  is  a  singularity  point 
of  zi  ),  and  around  each  new  singularity  C«j(^)  of  z^ .  (Note:  this  is  not  a  true  singularity 
of  z^,  as  there  is  smoothing  over  an  inner  scale  as  shown  in  section  7.)  The  scale  of  the  the 
inner  region  near  Co(f)  is  proportional  to  ,  whereas  the  difference  of  speed  of  the 

singularities  between  Co(0  C»i  is  order  unity.  This  means  that  for  t  =  0{B^I^)  , 

the  location  of  is  within  the  inner  scale  around  Co(0  aJid  the  two  term  (in  B  ) 
outer  expansion  that  one  can  expect  to  match  with  the  two  term  asymptotic  behavior  of 
the  solution  to  the  leading  order  inner  equation  for  large  values  of  inner  spatial  variable  is 
(from  Van  Dyke’s  matching  principle) 

*  ~  (Co(t).<)  (C  -  Cci(())  +  B  [  x,(()  (C  -  &(<))■’''  +  Mt)  (C  -  C„(<))‘“''] 

(4.10) 

The  matching  question  will  be  taken  up  in  section  6.  However,  if  indeed  a  successful 
inner  outer  matching  is  possible,  the  evidence  presented  in  this  section  suggests  that  new 
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singularities  for  zi  are  born  and  these  will  move  towards  the  physical  domain.  For  the 
case  of  SaiFman  solution,  each  such  newly  created  moves  faster  than  a  zero  since  we 
showed  in  section  3  that  the  speed  of  Co(0  smaller  than  for  this  class  of 

solutions.  We  emphasize  that  these  singularities  like  the  singularities  and  zeros  of  zo(^  are 
only  present  in  terms  of  the  outer  asymptotic  solution  (4.1)  and  that  they  are  smoothed 
out  in  an  inner  scale  as  will  be  seen  in  section  7.  In  section  7,  we  will  find  that  the  size  of 
each  such  inner  region  is  where  the  deviation  from  zqc  . 

Note  that,  if  we  were  to  look  at  the  evolution  equation  for  Z2  ,  it  would  contain 
singular  forcing  at  the  singularities  of  zq  and  the  zeros  of  •  In  addition,  it  will  have 
singular  forcing  at  the  singularities  of  zi  which  were  born  initially  from  a  zero  of  zq  ^  that 
moved  away  later.  Since  left  hand  side  of  the  differential  operator  in  (4.2)  is  the  same  for 
Z2  ,  Z3  ,  etc.,  as  it  is  for  zi  ,  there  are  no  new  singular  points  for  22  (or  ^3  ,  etc)  other 
than  ones  for  zi  . 

Note  that  the  singularity  structure  for  zi  and  22  and  other  higher  order  terms  that 
we  have  sketched  out  are  dependent  on  the  initial  condition  being  independent  of  B  .  If 
we  have  initial  conditions  that  allow  appropriate  terms  in  zi  ,  one  may  not  have  any  birth 
of  a  singularity  of  zi  from  an  initial  zero  of  Z(^  .  This  effect  will  show  up  at  22  or  higher 
order  depending  precisely  on  the  initial  conditions. 

Aside  from  each  point  Cs(0»  Co(0  and  corresponding  Csi(0>  where  the  perturbation 
expansion  (1.11)  fails  because  the  higher  order  terms  become  progressively  more  singular, 
(1.11)  is  also  invalid  at  any  singular  point  Cp(0  created  by  surface  tension  effects  (as  will 
be  seen  in  sections  5  and  6)  on  an  branch  point  or  zero  of  Zq(^,  where 

~  i„(()  (C  -  Cp(0)~''’  +  ■■  (4.11) 

as  C  — +  Cp(0  provided  where  Ao(0)  =  0  and  .<4o(t)  scales  with  some  some  positive  power  of 
B.  Indeed,  for  any  B,  which  need  not  even  be  small,  substitution  of  (4.11)  into  the  complete 
equation  (2.17)  (actually  its  ^  derivative)  shows  that  such  a  singularity  is  consistent  with 
the  equation  with 

i,  =  -«l(Cp(<).<)  -  (4.12) 

Because  of  the  scaling  of  Aq  with  some  positive  power  of  B,  the  singular  behavior  (4.11) 
is  not  present  in  2o<;;  further  since  Ao(0)  =  0,  there  is  no  such  singularity  at  initial  time. 
However,  at  later  times,  for  <  <<  B~^,  the  location  of  singular  point  (p{t)  that  initially 
coincides  with  some  C»(0)  (where  2^(^,0)  is  singular  as  in  (3.31)  with  0  <  <  2)  is  within 

an  inner  neighborhood  of  Ca(0  where  the  perturbation  expansion  (1.11)  is  invalid  in  any 
ca.se  because  of  the  nonuniformity  pointed  out  earlier.  Thus  there  is  no  need  for  a  separate 
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inner  region  around  such  Cp(0-  However,  in  section  (6),  we  shall  see  that  a  (^p{t)  which 
initially  coincides  with  with  some  zero  Co(0)  moves  away  from  Co{t)i  hut  for  t  =  0(1) 
remains  within  a  neighborhood  of  Csi(0-  However,  at  this  time,  the  the  inner 

region  around  C«i(^)  which  scales  lik  does  not  include  any  of  Cp(0*  A  separate 

inner  region  scaling  as  around  each  such  Cp(0  becomes  necessary  to  analyze. 

5.  Local  Nonlinear  Equations  near  singular  point  of  2:0 ^ 

Here  we  consider  the  equation  near  a  point  where  zo(;  is  singular  generally  like 

(C  -  Ut)r^  for  0  <  ^  <  2  . 

For  ^  >  2  ,  the  curvature  term  in  the  evolution  equation  is  less  singular  than  than 
other  terms  of  the  equation  and  there  is  no  need  for  an  inner  zone.  For  ^  <  0,  the  inner 
equations  becomes  a  little  more  complicated  and  is  not  considered  here.  The  equations 
derived  in  this  section  will  also  be  restricted  by  the  assumption  that  ICsit)]  —  1  is  not 
very  small.  Each  eventually  must  approach  |C|  =  1  as  shown  in  section  3  and  so 
our  assumption  implies  that  we  will  not  be  looking  at  the  final  stages  of  evolution  when 
|^,|  —  1  <<  1.  In  section  8,  however,  we  derive  equations  for  the  final  stages  of  evolution 
for  the  special  case  |  <  <  1  . 

Before  introducing  inner  variables,  it  is  convenient  to  take  the  derivative  of  (2.17)  with 
respect  to  Cj  to  come  up  with  an  evolution  equation  for  Z(.  Inorder  that  the  curvature 
term  comes  at  the  same  order  as  and  that  the  solution  matches  with  the  outer  solution 
for  which  zo(;  ~  jBo(t)(C  —  C«(t))~^?  it  becomes  necessary*  to  introduce  iimer  variables; 

C  -  C.(()  =  ^  X  (5-1) 

CM  G(x,  r)  (5.2) 

r  =  f  di  Cl'\t)  (C.  (t),  t)  (5.3) 

Jo 

where 

Co(t)  =  Co(0)  efo  (5.4) 

where  Co(0)  depends  on  the  initial  condition  Bo(0).  Then  the  leading  order  equation 
becomes: 

=  -2  ajj  (5.5) 

2 

In  deriving  the  above,  we  assumed  lC,(f)|  - 1  >>  ,  or  otherwise  the  simplification 

of  each  of  qi  through  q-j  by  ,..970  would  be  incorrect.  This  would  make  the  derivation 

*  The  term  C'o(t)  appearing  in  the  scaling  is  only  meant  to  put  the  inner  equation  in 
the  simplest  possible  form 
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of  (5.5)  invalid.  When  the  singularity  is  so  close  to  j^|  =  1 ,  different  scalings  are  necessary 
for  the  dependent  and  independent  variables  and  instead  of  a  local  partial  differential 
equation  (5.5)  for  the  inner  equation,  one  obtains  a  parameter  free  integro-differential 
equation.  We  consider  that  in  section  8  for  the  special  case  \  <  ji  <  \.  Equation  (5.5) 
is  the  Harry-Dym  equation  as  pointed  out  by  Kadanoff  (private  communiation)  which  can 
be  related  to  the  modified  KdV  equation  (Kawamoto,  1985)  and  is  completely  integrable. 
However,  all  the  inverse  scattering  theory  results  are  based  on  the  real  domain  where 
certain  decay  conditions  hold  at  ±oo  .  These  do  not  appear  to  carry  over  to  the  complex 
plane  and  we  could  find  no  direct  use  of  the  integrability  of  (5.5).  However,  there  exists 
similarity  solution  that  satisfy  both  initial  and  matching  conditions  at  infinity  for  certain 
ranges  of  arguments  of  Arg 

Consider  the  initial  value  problem,  with 

G(x,0)  =  i  (5.6) 

with  boundary  conditions 

fi(X.r)  ~  ~  (5.7) 

for  X  ill  some  sector  of  the  complex  x  plane  that  will  be  determined  shortly.  Equation 
(5.7)  corresponds  to  assuming  that  at  initial  time,  z,^  is  the  same  as  •  The  boundeiry 
condition  (5.7)  corresponds  to  the  leading  order  inner-outer  matching,  i.e.  as  we  move  out 
of  the  inner  zone  the  solution  is  assumed  to  match  to  zo(^  . 

It  is  convenient  to  transform  the  Harry  Dym  equation  (5.5) 

H  =  (5  8) 

Then  the  equation  for  H  is: 

Hr  =  (5.9) 

This  has  a  similarity  solution  that  satisfies  both  initial  and  boundary  conditions  (5.6)  and 
(5.7)  of  the  form 

H  (5-10) 

On  substituting  this  into  (5.9),  we  obtain 

1  ^  97  97  W 

g«^-2)(;9-4)  F  +  g(^-2)(2W"-84;?+80),  -  2fri‘  F"+j  (/?  -  2f  F'"  = 

(5.11) 

where 

^  ~  ^3-3^/2  (5.12) 
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Note  that  in  the  special  case  of  ^  =  1  ,  (5.11)  can  be  integrated  once  to  obtain 

3  _  27  9  27  3  1  1  1  /w  H  f%\ 

-r)F  -  —T]^  F'  -  —Ti^  F"  =  -[1  -  (5.13) 

where  we  used  the  botindary  condition  F  — »  1  as  t/  — >  0  . 

We  now  return  to  the  general  case  of  arbitrary  in  the  interval  (0,  2).  We  now  seek  the 
asymptotic  solution  to  all  powers  of  7/  in  order  to  use  Borrel  summation  to  find  the  leading 
order  transcendental  terms,  that  are  beyond  all  orders  of  the  asymptotic  expansion.  This 
will  give  information  on  what  sectors  of  the  complex  7/  plane  can  the  boundary  condition 
F{ri)  ~  1  be  satisfied  for  small  tj  corresponding  to  boundary  condition  (5.7)  for  some 
corresponding  sector  in  the  x  plane-  It  is  convenient  to  substitute 

F  =  1  +  Q  (5.14) 


into  (5.11)  to  obtain 

-E  (-1)”  5("+1)("+2)<3“0'  +  +  |(;8-2)(2W^-84;9+80).,Q' 

n=l 

0"  +  f  (^  -  2)’  v’  Q'"  =  O'  (5.15) 

It  is  clear  that  the  complete  asymptotic  expression  for  small  •q  is  of  the  form: 

OO 

Q  ~  am  (5.16) 

m=l 

On  substituting  this  into  (5.15),  we  obtain  a  complicated  recxirrence  relation.  However, 
in  order  to  find  the  leading  order  transcendental  terms  for  small  q  ,  we  only  consider  the 
form  of  the  recurrence  relation  for  large  tti  .  We  will  establish  aposteori  that  the  nonlinear 
terms  will  not  contribute  to  the  leading  order  in  the  ratio.  Consider  any  nonlinear  term 
in  the  summation  formulae  in  (5.15),  say  QQ'  .  The  coeflBcient  of  q”^~^  will  involve 
5  °m-i  •  From  the  final  asymptotic  expression  for  am  in  (5.19)  that  is 

obtained  by  balancing  the  linear  terms,  it  is  not  difficult  to  see  that  the  the  contribution 
from  the  previous  summation  in  equation  (5.15)  is  relatively  small  for  large  m  .  Similarly, 
we  can  argue  for  other  nonlinear  terms.  More  detailed  arguments  have  been  made  by 
Combescot  et  al  (1988)  along  these  lines  :br  a  nonlinear  ordinary  differential  equation 
arising  in  steady  state  selection.  While,  this  does  not  rule  out  other  forms  of  balance 
where  nonlinear  terms  can  be  important,  we  expect  that  the  balance  obteiined  through 
linear  equations  is  appropriate  since  for  large  \q\,  transcendental  terms  zure  small  in  certeiin 
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sectors  of  the  complex  plane  and  should  therefore  appear  as  solution  to  a  homogeneous  part 
of  equation  obtained  by  linearizing  about  F  =  1.  Kruskal  &  Segur  (1986)  demonstrated 
this  to  be  the  case  for  a  model  nonlinear  equation.  ^Prom  ignoring  the  nonlinear  terms  in 
(5.15),  we  obtain  the  recurrence  relation: 

-^  =  X  (^  -  2)’  (m  -  if  [I  +  0(i)l  (5.17) 

Ofn—l  o 

This  implies  that 

27 

~  l(m-l)!f  (5.18) 

where  b  depends  on  the  initial  few  iterates  where  nonlinearity  in  (5.15)  plays  a  role.  On 
using  Stirling’s  approximation,  (5.17)  simplifies  to: 


w27  .a”  (2m)! 


(5.19) 


where  b  =  y/ir  b  .  We  now  carry  out  the  procedure  of  Borrel  resummation,  as  previously 
done  for  the  steady  state  fingering  problem  by  Combescot  et  al.  Define  R(s)  by 


Ris) 


OO 


(5.20) 


Then  the  asymptotic  behavior  (5.19)  implies  that  the  nearest  singularity  of  R{s)  to  the 
origin  in  the  complex  s  plane  is  at  s  =  —sq  ,  where  sq  =  and  the  leading 

order  asymptotic  behavior  near  s  =  —sq  is  given  by 


R{s)  ~  h  r(-l/2)  (s  + 


(5.21) 


The  functions  R{s)  and  F{ri)  are  related  through  the  transform: 


Alternately,  we  can  write 


ds  e  *  R{s'^  rj) 


(5.22) 


1  /*oo 

Qin)  =  ^  /  dsie-^  R{sl)  (5.23) 

yV  Jo 

Using  Watson’s  Lemma,  it  clearly  follows  that  the  asymptotic  expansion  of  Q  is  as  given 
by  (5.16)  for  r)  real  positive  and  small.  This  is  also  true  for  Arg  77  in  (— tt,  tt)  since  the 
contour  along  the  real  axis  is  not  affected  by  the  presence  of  singularity  of  R{s)  .  Notice 
that  instead  of  (5.22)  or  (5.23),  which  are  equivalent,  we  could  have  equally  used 


(5.24) 
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for  any  real  and  positive  constant  c  that  is  independent  of  77  and  obtain  the  same 
asymptotic  expansion  (5.16)  for  Arg  77  in  (—ir,  x)  .  It  is  clear  that  (5.23)  and  (5.24) 
differ  by  transcendentally  small  terms  in  77  that  axe  beyond  all  orders  of  the  asymptotic 
expansion  (5.16).  However,  we  shall  see  momentarily  that  (5.24)  does  not  give  the  asymp¬ 
totic  expansion  (5.16)  when  Arg  Tf  is  extended  beyond  x  which  is  necessary  for  obtaining 
matching  to  the  outer  asymptotic  expansion  for  a  range  of  P  values. 

Returning  to  (5.22),  as  we  continuously  change  argument  of  77  from  0  to  tt  ,  the 
stationary  phase  path  encounters  the  singularity  at  s  =  Sc  =  i  on  the  real  positive 
axis  when  Arg  77  =  tt  ,  In  conformity  of  emalytic  continuation  from  smaller  values  of 
Arg  77  ,  it  follows  that  the  path  of  integration  should  be  deformed  such  that  it  goes  below 
Sc  on  the  real  axis.  For  2  x  >  Arg  77  >  tt  ,  this  deformation  is  equivalent  to  Fig. 

6,  which  is  the  required  stationary  phase  contour.  The  contribution  from  the  trip  around 
77  =  Sc  gives  an  additional  contribution  of 

e*  2^/2  b  X  77^/^  So  e"‘V^  (5.25) 

However,  this  is  transcendentally  small  for  Arg  77  in  the  interval  (tt,  2  x)  and  hence 
the  asymptotic  expansion  (5.16)  holds  for  this  range  of  Arg  77  .  Once  again,  if  we  repeat 
the  analysis  for  Arg  77  continuously  changing  from  0  to  —x  and  going  beyond  —x  ,  we 
find  that  there  is  an  additional  contribution  from  s  =  ~Sc  owing  to  the  use  of  a  contour 
shown  in  figure  7  that  is 

—  e~‘  2?!^  b  X  77^/^  e*V^  (5.26) 

However,  this  is  transcendentally  small  for  Arg  77  in  the  interval  (—277,  — tt)  ,  Thus  over 
Arg  Tf  in  the  range  (—277,  277)  ,  the  asymptotic  series  for  Q  in  (5.16)  hold  for  small 
77  implying  a  matching  is  possible  in  that  sector  of  the  complex  77  pleme  to  the  appropriate 
zero  surface  tension  solution  eis  (  —7  Ca(0  •  However  for  Arg  77  in  (2  77,  3  77)  and 
(—3  77,  —2  77)  ,  (5.16)  does  not  hold  as  (5.25)  or  (5.26)  become  transcendentally  large.  At 
this  stage,  we  can  see  the  difference  between  (5.24)  and  (5.22)  (or  (5.23)).  If  we  were  to 
use  (5.24),  the  asymptotics  of  Q  when  Arg  r]  is  in  (77,  277)  or  (-277,-77)  will  be 
dominated  by  the  contribution  from  c  which  would  then  dominate  (5.16).  If  c  were 
complex  with  Re  c  >  0  ,  it  would  still  not  be  possible  to  ensure  the  validity  of  the 
asymptotic  series  (5.16)  beyond  an  interval  of  2  77  for  Arg  77  .  We  now  note  that  there  is 
nothing  special  about  the  interval  (—277,  277)  for  Arg  77  where  the  asymptotic  series  (5.16) 
is  valid.  Through  appropriate  rotation  of  contours  in  the  s—  plane  in  the  expression  (5.22), 
(5.16)  can  be  made  valid  for  any  interval  of  the  form  (—277  -|-  2  fc  77,  2  fc  77  -|-  2  77)  for  any 
integer  k. 


Now,  we  examine  the  implications  of  being  able  to  carry  out  an  inner-outer  match  for 
a  range  of  Arg  rj  not  exceeding  47r.  Clearly,  the  matching  fails  because  the  transcendental 
term  (5.25)  or  (5.26)  becomes  large.  These  transcendental  behavior,  when  rewritten  in 
terms  of  outer  variable  correspond  to  transcendental  corrections  to  the  regular  perturbation 
expansion  (1.11)  of  the  form 

(5. 


where  in  (5.27),  we  have  suppressed  a  prefactor  that  is  algebraic  in  B  and  C  —  Ca(0-  This 

2 

form  of  transcendental  correction  to  ( 1 . 1 1 )  is  only  valid  when  B <<  —  Cs(t)|<<l. 

To  find  the  form  of  the  transcendental  terms  for  arbitrary  C>  one  needs  to  look  at  the 
associated  homogenous  equation  found  by  linearizing  (2.17)  about  ^o  and  look  for  WKB 
type  solution  to  this  equation.  Each  WKB  solution  obtained  this  way  are  expected  to  have 
the  feature  that  they  exchange  dominancy  across  lines  in  the  complex  ^  plane  called  Stokes 
lines.  The  structure  of  global  Stokes  lines  fo  the  steady  finger  problem  has  been  analyzed 
before  (Tanveer  (1987a)).  However,  in  the  time  dependent  problem,  the  equations  appe^lr 
to  be  quite  complicated  and  we  do  not  address  this  problem  here. 

We  now  present  just  heuristic  arguments  to  conclude  that  appropriate  matching  can 
be  accomplished  atleast  for  some  range  of  /3.  These  have  to  be  checked  in  the  future  by 
complete  determination  of  global  Stokes  lines  at  each  instant  of  time.  ^From  prior  expe¬ 
rience  with  steady  problems,  it  appears  likely,  that  in  |^|  =  1,  the  dominant  contribution 
towards  the  transcendental  correction  is  from  the  nearest  singularity.  It  is  clear  from  (5.27) 
that  the  local  structure  of  the  Stokes  lines  near  C  =  C»(^)  is  given  by  setting  the  real  part 
of  the  exponent  in  (5.27)  to  zero.  Now,  from  the  relation  (5.1),  (5.3)  and  (5.12),  it  follows 
that 


Arg 


(C.(0  -  C) 
Ut) 


Arg  Csjt) 
(3  -  3^/2) 


-  Arg  Co{t)  - 


1 

3  -  3/3/2 


Arg  rj 


(5.28) 


Now,  in  the  special  case  of  a  /9  =  1,  i.e.  a  simple  pole  for  zo(^,  it  is  clear  that  regardless 
of  the  value  of  tt  —  -  Arg  C„«),  ein  integer  k  can  be  chosen  so  that  Arg  rj  is  in 

(— 27r  -f  2k7T,  2n  4-  2kn)  when  Arg  is  in  (— 7r/2,  7r/2)  showing  that  a  matching 

is  possible  to  the  outer  solution  in  eill  directions  from  (g{t)  towards  the  half  plane  that 
contains  the  physical  domain  (see  Fig.  8).  We  now  argue  that  such  a  matching  imposed 
for  the  nearest  singularity  to  the  physical  domain  |C|  <  1,  assumed  to  be  at  0(1)  distance, 
is  sufficient  to  ensure  the  validity  of  (1.11)  in  the  physical  domain.  Consider  a  circle  around 
the  origin  in  the  (  plane  so  that  its  radius  R  is  slightly  less  (independent  of  B)  than 
the  distance  of  the  nearest  singularity  5  from  the  origin  (See  Fig.  8).  We  know  from 
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our  discussions  in  the  context  of  zero  surface  tension  solution,  that  the  value  of  z  on  this 
large  circle  uniquely  determines  z  on  |C|  =  1  in  a  well  posed  manner.  Now  if  any 
part  of  this  large  circle  contained  exponentially  large  terms  in  B  that  can  be  expected 
if  an  inner-outer  matching  failed,  then  one  or  more  of  the  coefficients  of  a  power  series 
representation  for  z  -f-  ^In  C  for  the  channel  case  {  z  —  for  the  radial  geometry) 
that  is  convergent  up  to  this  circle  of  radius  R  must  contain  transcendentally  large  terms 
in  B  .  Now  since  the  distance  of  the  singularity  from  the  unit  circle  is  assumed  to  be  order 
unity  (i.e.  independent  of  6  ),  it  follows  that  z  on  |C|  =  1  will  contain  transcendentally 
large  terms  in  B  that  is  incompatible  with  the  assumption  that  (4.1)  actually  holds  in  the 
physical  domain. 

The  argument  just  presented  above  suggests  that  for  a  simple  pole  of  zo(^  ot  j3  values 
close  to  unity,  we  are  assured  of  a  successful  inner-outer  matching.  However,  this  condition 
may  just  be  sufficient,  not  necessary.  Certainly,  the  argument  above  will  not  hold  for  any 
C*(t)  which  is  not  the  nearest  singularity  from  the  physical  domain.  Further,  the  global 
structure  of  the  Stokes  line  is  likely  to  change  due  to  effect  of  other  singular  points  that 
move  out  of  the  inner  region,  if  these  are  located  close  to  a  Stokes  line  emanating  from 
Cj(<).  This  point  is  not  addressed  in  this  paper,  but  needs  to  be  resolved  in  the  future 
to  ensure  that  the  inner-outer  structure  of  the  solution  is  indeed  valid  for  any  ^  in  the 
interval  (0, 2).  Our  argument  only  assures  that  for  /?  =  1  or  sufficiently  close  to  it. 

Now  consider  possible  singularities  of  the  actual  Z(^.  It  is  clear  that  the  function  F 
appearing  in  the  similarity  solution  (5.10)  that  satisfies  (5.11)  admit  singularities  at  one 
or  more  tjo  such  that 

F{v)  ~  ii  in-Vof"  (5.29) 

as  r;  7/0.  Rewritten  in  terms  of  the  outer  variable,  this  implies  that 

~  (^)’^ir'(3-3/)/2r'‘/’  IC-CpW]'"'’  (5.30) 


for 

lC-Cp(t)l  «  IlCoCf))-  [— 

where 

Bt 

Cpit)  =  Cs(t)-{Co{tr^  — 

v> 

We  now  present  a  likely  explanation  to  Dai  et  eiI’s  (1991)  numerical  computation  that 
appeared  to  suggest  that  an  initial  pole  of  ^<(C>0)  remains  a  pole  at  later  times  even  for 
nonzero  B  that  is  small.  In  order  to  determine  the  nature  of  the  singularity  Cp(0  (5.30) 
on  |(^|  =  1  through  monitoring  of  the  power  series  coefficient  in  (1.6)  as  carried  out 
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by  Dad  et  al  (1991),  one  must  consider  the  asymptotic  behavior  of  A;„  for  very  lairge  n 
because  of  the  highly  localized  nature  of  the  singularity  as  reflected  in  (5.31).  Further,  the 
calculations  must  be  accurate  enough  to  distinguish  the  the  small  0(B^/®)  (for  /3  =  1) 
coefficient  of  the  singular  term  in  (5.30)  from  the  behavior 

Bo 

that  is  valid  for  «  |C  —  C«(0l  ^  given  that  ^*(0  and  Cp(t)  are  quite  close  to 

ea.ch  other  (within  as  is  seen  in  (5.32)).  We  believe  that  the  behavior  (5.30)  will  be 
reflected  in  in  fcn  for  very  large  n,  far  beyond  what  has  been  calculated  by  Dai  et  al  (1991). 
In  some  sense,  in  the  Dai  et  al  (1991)  calculation,  the  singular  behavior  (5.33)  masks  the 
actual  singularity  because  sufficiently  small  spatial  scales  were  not  resolved. 

However,  due  to  the  stiffness  of  the  system  of  ODE  for  small  B  and  larger  size  of 
truncation  in  the  Dai  et  al  (1991)  formulation,  such  a  direct  computation  appears  to  be 
impractical  at  this  stage.  This  The  number  of  kn  numerically  computed  by  Dai  et  al’s 
(1991)  does  not  appear  to  be  large  enough  and  accurate  enough  to  detect  the  singular 
behavior  (5.30).  Instead,  their  calculations  only  reflected  the  behavior:  that  is  valid 
6.  Local  Nonlinear  Equations  near  a  zero  of  zo( 

As  pointed  out  earlier,  the  perturbation  expansion  (4.1)  is  invalid  near  the  zero  (o{t) 
of  zo^  as  well.  We  now  present  the  analysis  for  an  inner-region  around  (o{t)  which  scales 
as  where  (4.1)  breaks  down. 

We  introduce  innter  variables 


where 


C  -  Co(t)  =  6"/"  h  i 

0-2/7  f*  g2o^(Co  jt'),  t') 

Jo  ^2  (t') 

=  H"/"  k2  G  (e,  r)  j*  dt'  92(Co(<')^  <') 


^  _ gy(Co(<)>0 _ 

92of^(Co(t),<)  2^0 Jf  (Co(0.0 

92oJ^^(Co(t),t) 

Then  equation  (2.15)  to  the  leading  order  in  B^^^  reduces  to 


ko  = 


=  -G(  +  ^  -  2 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 
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In  order  that  this  matches  with  as  C  Co(t)  ?  we  must  require  that 


G^(e,r)  ~  i  (6.7) 

as  ^  »  oo  along  sectors  in  the  complex  ^  plane  for  any  fixed  scaled  time  t  in  directions 

towards  the  physical  domain.  The  initial  conditions  imply  that 

G^(^,0)  =  i  (6.8) 


For  convenience,  we  define 

H  =  (6.9) 

Then  on  taking  the  derivative  of  (6.6)  with  respect  to  ^  ,  we  find  that  the  evolution 
equation  for  H  is  given  by 

-2  -  2  ^  =  1  -  2Hut  (6.10) 

The  boundary  and  initial  conditions  (6.7)  and  (6.8)  translate  to: 

~  r"'"  (6-11) 


for  1^1  >>  1  and  initially 

^(^,0)  =  r'/"  (6.12) 

6a.  Early  times  t  « 

When  t  «  ,  it  is  clear  from  the  definition  of  t  in  (6.2)  that  r  <<  1  .  In 

this  case  there  is  an  asymptotic  similarity  solution  that  is  relevant  to  boundary  and  initial 
conditions  and  is  given  by 

(6.13) 

where  F{ri)  (r/  =  y,;,  in  this  ceise)  satisfies  (5.11)  with  /3  =  — 1  and  can  be  made 

to  satisfy  boundary  condition 

lim  F{r])  =  1  (6.14) 

for  Arg  t/  in  (2fc  tt  —  27r,  2kTT  +  2  tt)  for  any  specific  choice  of  integer  k,  as  shown  in 

section  5.  Using  the  relations  (6.1),  (6.2)  and  the  relation  of  7/  to  r  and  ($  —  r,  it  is  cleeir 

that  this  implies  that  the  inner-outer  matching  condition  appearing  as  boundary  condition 
(6.11)  will  be  valid  for  a  Stt/Q  range  of  Arg  ^  and  hence  of  Arg  (1  —  C/Co(0). 
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We  now  determine  the  singularity  locations  of  F(r/)  for  a  solution  that  satisfies  (6.14) 
for  Arg  rj  in  the  intervEil  {—2ir,  27r).  Through  appropriate  rotational  invEiriances  of  the 
equation  (6.6)  and  initial  conditions,  the  form  of  the  singularities  for  solutions  that  satisfy 
(6.14)  for  Arg  rj  in  {—2ir  +  2kir,  27r  +  2A:7r)  for  nonzero  integer  k  C2m  then  be  deduced.  The 
proper  choice  of  integer  k  depends  on  Arg  ki  and  Arg  ^2  as  determined  in  (6.4)  and  (6.5) 
as  they  relate  the  Arg  (1  —  C/Co(0)  ■^‘''9  ^  hence  to  Arg  tj. 

First,  for  computational  purpose,  we  found  it  convenient  to  introduce  the  transforma¬ 
tion 

T]  =  (6.15) 

with  principal  choice  of  Arg  tj  and 

F(ri) = (sie) 
Then  the  differential  equation  satisfied  by  F  is 

-5^  [^  +  2  ^  ^']  = 

The  boundary  condition  (6.14),  imply  that 

Hv)  ~  (6-18) 


Through  consistent  dominant  balance  argument  on  (6.17),  one  obtains  the  following  pos¬ 
sible  eisymptotic  series  at  00: 


225  X  115  _ 19/2 


(6.19) 


The  form  of  the  transcendental  correction  has  already  been  determined  in  section  (5)  for 
F{i])  as  T)  — »  0.  Using  that  result,  we  can  assure  ourself  that  there  is  not  transcendental 
term  (atleast  to  the  leading  order)  when  Arg  fj  is  in  the  interval  (— 27r/9,27r/9).  At  the 
anti-Stokes  lines  27r/9  and  — 27r/9,  transcendental  terms  are  born  so  that  for  Arg  rj  in  the 
interval  (27r/9,67r/9), 


^-1/2  _  ^  ^  ^  ^  (J,  fi-13/8^W32/729ii^/* 


(6.20) 


where 

Cx  =  (6.21) 


and  for  Arg  fj  in  (— 67r/9,  — 27r/9) 

F  ^  fj-^/^  -  +  225j^  --19/2  ^  ^  (g  32) 

8  128 
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where 


Ca  =  5"^/^  (6.23) 

Notice  that  the  transcendental  terms  in  (6.20)  and  (6.22)  are  subdominant  to  the  algebraic 
terms  of  the  asymptotic  expansion  for  Arg  fj  in  (2?r/9,  4t/9)  and  (— 47r/9,  —2  ir/9)  respec¬ 
tively,  but  become  the  dominant  contribution  for  Arg  fj  when  the  Stokes  line  Arg  fj  =  47r/9 
or  Arg  fj  =  — 47r/9  is  crossed.  Thus  (6.19)  holds  for  Arg  fj  in  (— 47r/9,  47r/9)  since  the 
transcendental  terms  in  (6.20)  and  (6.22)  axe  subdominant. 

We  now  describe  a  numerical  procedure  to  calculate  such  a  solution  F  to  (6.17)  in  the 
complex  fj  plane  and  to  find  the  location  of  singularities  fjo,  where  F  is  singvdax.  Earlier, 
Constantin  &  KadanofF  (1991)  and  Howison  (1991)  realized  that  a  two-thirds  singularity 
of  F  was  consistent  with  (6.17).  Indeed,  an  expansion 

F  ^  [fj-fjo]^^^  (6.24) 

i=o 

is  consistent  with  (6.17).  We  therefore  proceeded  to  verify  that  the  form  (6.24)  is  correct, 
determine  the  location  of  each  such  singular  point  and  find  out  if  there  are  other  forms 
of  singularities  admitted  by  (6.17). 

As  suggested  originally,  by  Kruskal  &  Segur  (1986)  for  another  third  order  nonlinear 
ODE,  a  convenient  numericed  method  to  calculate  such  solutions  will  be  to  march  in  along 
the  positive  real  axis  from  a  large  distance  by  using  the  asymptotic  behavior  (6.19)  for  F, 
F'  £ind  F".  This  ensures  that  we  have  effectively  set  the  coefficients  of  the  two  possible 
exponential  terms  to  zero  for  Arg  fj  in  the  interval  (— 27r/9,2  7r/9),  since  otherwise  on 
the  real  positive  axis  (a  Stokes  line),  each  of  the  two  exponential  exponential  contributions 
would  have  dominated  some  algebraic  terms  that  are  included  in  the  three  term  asymptotic 
expansion  (6.19).  The  third  degree  of  freedom  for  the  third  order  ODE  (6.17)  has  been 
utilized  in  demanding  (6.18)  rather  than  the  more  general  condition  F  ~  constant  fj~^^^, 
which  for  arbitrary  constant  is  a  possible  leading  order  asymptotic  behavior  for  a  solution 
to  (6.17).  We  did  not  use  a  starting  value  of  fj  exceeding  15  to  avoid  the  expected  stiffiiess 
of  the  differential  equation  for  large  fj.  However,  for  a  starting  fj  exceeding  10,  we  verified 
that  the  computed  solution  F  was  independent  of  starting  to  10  digits.  A  Runge-Kutte 
solver  that  automatically  adapts  the  step  size  to  control  errors  was  to  calculate  numerical 
solution  for  any  given  fj  «  w  the  real  axis  for  fj  in  the  interval  (—5,15).  This  numerical 
integration  showed  no  sign  of  singularities  on  the  real  axis  on  this  interval  and  F  was 
always  real  upto  machine  precision.  We  then  used  a  procedure  quite  similar  to  one  used 
earlier  for  computing  singularities  in  a  different  context  (Tanveer,  1991)  We  took  round 
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trip  paths  over  iinti-clockwise  closed  contours  C  in  the  fj  plane  and  numerically  evaluated 
the  values  of 


(6.25) 


1 

Am 


(6.26) 


and  found  it  was  zero  to  numerical  acciuracy  most  of  the  times,  suggesting  that  there  are 
no  singularities  inside  such  contours  C.  Further,  it  was  checked  that  F  returned  to  the 
same  value.  When  this  contour  integration  answers  were  nonzero  for  any  contour,  we  also 
fovmd  that  F  during  those  times  did  not  return  to  the  same  value.  So  we  suspected  that 
there  was  one  or  more  singularities  inside  our  contour.  In  that  case,  we  went  around  that 
contour  C  anti-clockwise  three  times  and  evaluated  each  of  the  integrals  (6.25)  and  (6.26). 
Then,  we  either  found  that  (6.25)  gave  us  a  value  of  1  (within  numerical  accuracy)  and 
during  those  times  F  returned  to  the  same  value  after  three  rotations;  or  at  other  times  the 
value  of  (6.25)  was  not  unity  and  during  those  times  F  did  not  return  to  the  same  value 
even  after  three  rotations.  In  the  first  case,  we  concluded  that  there  was  one  singularity  of 
the  form  (6.24)  at  a  location  fjo  given  by  the  value  of  the  integral  in  (6.26)  (remembering, 
offcourse,  that  C  is  traversed  thrice).  Such  a  value  was  consistent  with  the  location  of  the 
actual  contour  C,  regardless  of  its  size,  which  we  varied.  In  the  second  case,  we  concluded 
that  the  contour  included  mtiltiple  singularity  (which  need  not  be  on  the  same  Riemann 
sheet).  In  such  cases,  we  decreased  the  size  of  the  contour  C,  imtil  we  found  that  the  value 
of  (6.25)  on  a  thrice  traversed  contour  was  unity. 

With  our  numerical  search  for  ^  in  the  region —5  <  Re  ij  <  5,  —7.5  <  Imi]  <  7.5, 
we  were  only  able  to  find  three  pairs  of  singularities  located  at 


fjo  =  0.29349  ±  3.7586  i,  0.660998  ±  5.61650  i,  and..  0.92209  ±  6.9174  i  (6.27) 


It  was  difficult  to  get  reliable  values  of  the  integrals  in  (6.25)  and  (6.26)  for  larger  values 
of  fj  far  away  firom  the  real  axis  because  of  the  effect  of  exponentially  growing  terms.  This 
led  to  restrictions  on  the  zone  where  reliable  singularity  search  could  be  performed.  One 
curious  observation  about  the  roots  in  (6.27)  are  that  they  are  close  to  the  Stokes  line 
Arg  fj  =  ±4  7r/9  that  limit  the  validity  of  the  asymptotic  exp2insion  (6.19)  and  get  closer 
as  the  roots  are  further  out  from  the  origin.  Now,  a  divergent  asymptotic  series  for  F  valid 
in  some  sector  of  the  complex  plane  has  encoded  in  it  information  about  the  asymptotic 
behavior  for  other  sectors  where  the  zisymptotic  series  is  invalid  due  to  its  subdominance 
to  transcendental  terms  that  now  become  large.  However,  the  form  of  the  transcendental 
correction  when  Arg  fj  is  just  beyond  47r/9  is  given  the  one  in  (6.20).  Now  notice  that 
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the  singularities  in  (6.27)  axe  at  reasonably  large  distances  from  the  origin,  especially  the 
latter  ones.  We  argue  that  a  singularity  fax  out  from  the  origin  ought  to  be  located  close 
to  a  point  fjg,  where 


+  Cifj, 


(6.28) 


Offcourse,  the  asymptotic  series  (6.20)  is  itself  valid  in  the  immediate  vicinity  of  such 
a  point  fjg,  as  can  be  expected  since  F  satisfying  (6.17)  does  not  admit  a  simple  zero. 
However,  we  mahe  the  ansatz  that  the  modification  of  (6.20)  occurs  over  a  small  local  scale 
when  ijo  is  large  and  so  fjg  is  close  to  rjo.  This  is  admittedly  a  very  heuristic  argument; 
however,  we  now  see  the  ramifications  when  this  is  assiuned  correct.  Consider  the  zeros  of 
(6.28)  near  the  Stokes  line  Arg  fj  =  47r/9.  We  put 


fjg  =  r  +  ‘‘ 

and  carry  out  a  perturbation  for  large  r  to  find  that 


9v^r9/4  [8 


-  Inr  -  In  |Ci| 


1  ro  ,  3  tt//® 


Computing  (6.31)  for  fc  =  1,  2,  3  gives  us: 


r  =  3.6788,  5.6252,  and  6.961618 


(6.29) 


(6.30) 


(6.31) 


(6.32) 


which  compares  well  with  the  actual  values  of  |^ol  for  the  roots  given  by  (6.27): 


3.77004,  5.65526..,  and  6.97858  (6.33) 

The  agreement  between  (6.32)  and  (6.33)  gets  better  for  larger  singularity  distances  from 
the  origin,  as  can  be  expected  if  our  heuristic  argument  is  correct.  Now,  lets  compare 
the  value  of  e  for  k  =  1,2,3,  with  what  is  observed  to  be  the  for  deviation  of  Arg  fjo 
from  47r/9  for  the  roots  in  (6.27)  in  the  upper  half  plane.  We  cannot  directly  apply  the 
asymptotic  relation  (6.31)  accurately  since  \Ci  \  given  in  (6.21)  is  only  known  interms  of 
a  real  posive  number  b  (introduced  in  section  5)  that  we  have  never  computed.  Leaving 
out  the  ICi  I  term  in  (6.30)  would  have  been  allright  as  a  leading  order  approximation  if 
r  was  sufficiently  large  so  that  \ln  r\  »  1.  But  this  is  not  the  case.  We  circumvent 
this  difl&culty,  by  estimating  \C\  \  by  using  the  last  value  in  (6.33)  and  using  it  in  (6.30)  to 
determine  |Ci|.  This  computed  \C\  \  is  then  used  in  (6.30)  for  k  =  1,  2  to  give  us: 


c  =  0.095909,  0.05768 


(6.34) 


compared  to  0.096606  and  0.05738,  respectively  for  Arg  fjo  —  47r/9  for  the  first  two 
roots  quoted  in  (6.27)  (only  those  in  the  upper  half  plane).  Again,  with  the  asymptotic 
axgtunent  results  is  better  the  further  out  the  singularity  is.  Thus,  it  appears  that  our 
heuristic  argument  about  the  infinite  set  of  singularities  along  the  Stokes  line  4ir/9  may 
be  correct.  If  so,  their  radial  positions  approach  those  given  by  (6.31)  for  to  positive 
integral  values  of  k  in  (6.31)  (the  agreement  getting  better  for  larger  k),  whose  Arg  values 
approach  47r/9  from  above  with  the  deviation  e  approaching  that  predicted  by  (6.31).  We 
checked  the  k  =  4  case  with  direct  ntunerical  computation  which  yielded  a  singtdar  point 
at  ijo  =  (1.128,7.966),  very  close  to  what  is  predicted.  Clearly,  from  symmetry  of  the 
problem  about  the  positive  real  axis,  on  the  basis  of  our  heuristic  argiunent,  there  will 
be  an  infinite  number  of  roots  close  to  the  Stokes  line  Arg  fjo  =  — 47r/9,  with  complex 
conjugate  locations  to  the  ones  determined  above. 

Now  let’s  examine  the  consequence  of  singularities  in  the  similarity  solution.  Going 
back  to  the  relation  (6.13),  (6.15)  and  (6.16)  it  is  clear  that  for  r  <<  1, 

m,r)  ~  -r-’’!”  (6.35) 

for 

l«  -  «,(’•)!  «  (6.36) 

where 

4p(r)  =  r  -  (6.37) 

Thus  initally  each  singularity  coincides  with  C  =  0,  i.e.  Co(0).  However,  they  move 
away  from  it  for  r  >  0.  Now  for  any  r  >  0,  assuming  that  there  are  infinite  number 
of  singularities  that  are  arbitrarily  far  out  in  the  rj  plane  eis  indicated  above  that  will 
approach  the  Stokes  lines  that  limit  the  validity  of  the  eisymptotic  expansion  (6.19)  (and 
hence  inner-outer  matching).  These  singular  points  located  in  |^|  >>  1  are  expected  to 
alter  the  global  Stokes  line  structure  for  C  outside  the  inner  region  around  Co(0’ 
unclear,  what  this  affect  will  be.  However,  we  proceed  further  with  the  assumption  that 
that  the  global  Stokes  line  structure  is  such  that  local  matching  for  Arg  fj  (and  hence 
Arg  (  for  large  ()  in  the  interval  {~4kI9,  47r/9)  guarantees  that  transcendentally  large 
terms  in  B  are  avoided  on  1(|  =  1  thereby  ensuring  the  validity  of  (1.11). 

6b.  The  case  of  f  = 

In  this  case  r,  as  defined  in  (6.2),  is  order  unity.  The  similarity  solution  (6.13)  becomes 
invalid  and  one  has  to  look  at  other  ways  of  solving  the  P.D.E.  (6.10).  Our  primary  interest 
in  (6.10)  is  to  see  if  solutions  exist  that  satisfy  the  boundary  condition  (6.11)  over  some 
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sector  of  the  complex  ^  plane,  if  for  early  times  such  a  solution  can  be  matched  to  the 
early  time  similarity  solution  (6.13)  and  what  singularities  are  present  in  such  solutions. 

We  assume  that  the  only  singularities  of  (6.6)  are  those  created  at  r  =  0,  though 
located  at  a  different  position  that  is  determined  from  (4.12).  Note  that  near  a  singularity, 
the  ^  term  has  little  importance  and  therefore  the  equation  in  an  approximate  sense  reduces 
to  a  Harry-Dym  equation  (with  the  independent  variables  ^  —  r  and  r)  for  which  it  might 
be  expected  that  the  only  moving  singularities  are  of  -4/3  rd  type  for  G^.  With  this 
assumption,  we  expect  an  infinite  number  of  singularities  moving  away  from  ^  =  0.  Now 
consider  the  question  of  inner-outer  matching. 

What  is  necessary  is  to  ensme  that  boimdary  conditions  (6.11)  can  be  satisfied  for 
sufficiently  large  range  of  Arg  For  that  piurpose,  we  need  to  look  for  other  possible 
behavior  of  H  aside  from  terms  that  are  algebraic  in  ^  that  occur  in  systematic  domi¬ 
nant  balance  arguments  applied  to  (6.10)  and  (6.11).  This  can  be  found  by  linearizing 
(6.10)  about  the  asymptotic  behavior  (6.11)  and  looking  for  solutions  H  to  the  associated 
homogeneous  equation,  which  in  this  case  turns  out  to  be: 

+  (6.38) 

We  impose  the  condition  that 

HU,0)  =  0  (6.39) 

for  einy  ^  ^  0,  as  it  appears  appropriate  beyond  the  immediate  scale  where  the  similarity 
solution  of  section  6a  is  relevant.  Also,  as  ^  oo  for  certain  range  of  arguments  that  is 
\mspecified  at  this  time, 

0  (6.40) 

Such  a  solution  can  be  found  in  the  form 

H  =  e^  (6.41) 


The  leading  order  WKB  type  behavior  of  W  is  given  by  Wq  that  satisfies 


[Wor  +  W„(]  -W„l  =  0 


(6.42) 


Appropriate  solution  to  this  which  matches  with  the  form  of  the  exponential  in  (6.20)  and 
(6.22)  (after  appropriate  change  of  variables)  for  r  <<  1  is  is 


Wo  =  ± 


iVg2(4  -  rf* 

t>/29’/’ 


(6.43) 
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The  next  order  correction  to  W  that  matches  with  the  form  of  the  algebraic  prefactors  in 
(6.20)  and  (6.22)  (after  appropriate  change  of  variables)  is 


One  can  clearly  see  that  the  form  of  the  leading  order  transcendental  corrections 
agrees  with  with  the  form  for  r  <<  1,  as  reflected  in  (6.20)  and  (6.22). 

iProm  the  (6.41)  and  (6.43),  it  is  clear  that  with  a  choice  of  positive  sign,  is  tran- 
scendentally  small  for  large  ^  when  Arg  ^  is  in  the  interval  (0,  |7r)  and  with  the  minus  sign 
choice  ff  transcendentally  large  for  large  ^  in  this  sector  but  transcendentally  small  when 
Arg  ^  is  in  the  adjoining  interval  (— |7r,0).  Now,  for  a  general  partial  differential  equation 
with  third  order  spatial  derivative,  one  can  expect  to  be  able  to  specify  three  independent 
boundary  conditions.  For  (6.10),  one  degree  of  freedom  has  been  used  in  requiring  (6.11) 
rather  than  the  more  general  asymptotic  condition  H  ~  [constant  (^  — 'r)  +  T]~^^^. 
The  other  two  degrees  of  freedom  can  be  used  by  ensuring  that  the  the  coefficients  of 

(^Irjia/a ^  for  Arg  (  in  (— |7r,0]  and  and  coefficients  of 

for  Arg  ^  in  [0,  |7r)  are  all  zero.  For  r  <<  1,  we  have  seen  in  terms  of  compuation  of  a 
similarity  solution,  that  these  conditions  are  appropriate. 

This  would  mean  that  it  is  possible  to  satisfy  the  asymptotic  boundary  condition 
(6.11)  for  Arg  (  in  the  open  interval  (— ^,  x)*  There  is  nothing  special  about  this  range 
of  argument  of  Clearly,  instead  of  this  range,  we  can  choose  to  ensure  boundary  condition 
(6.11)  for  other  range  of  arguments  like  (0,  ^)  ),  etc.  The  appropriate  choice  will  depend 
on  97o(Co(^)iO»  relate  Arg  (C  —  Co(^))  to  ■^'"'9  ^  through  (6.1).  Certainly  |7r  is 

smaller  than  tt.  Thus  the  proposed  sufficient  condition  (section  5)  for  inner-outer  matching 
for  the  neeirest  singularity  is  not  applicable  here.  However,  if  we  are  to  assume  that  there 
are  an  infinite  set  of  —4/3  singularities  (as  our  arguments  indicate)  that  change  the  global 
Stokes  lines  in  such  a  way  that  the  matching  in  an  appropriate  Stt/Q  ensures  the  validity 
of  (1.11)  in  the  physical  domain,  then  we  have  a  consistency  in  inner-outer  matching.  We 
cannot  ascertain  this  at  this  stage  without  studying  the  global  Stokes  line  structure,  but 
proceed  further  on  the  assumption  that  this  is  the  case. 

Further  in  the  sectors  where  the  asymptotic  expansion  (6.11)  is  valid,  an  examination 
of  the  two  term  asymptotic  relation  for  G(  =  ,  obtmned  directly  from  (6.6)  and 

using  by  using  the  leading  order  relation  (6.7)  implies  that 


(6.44) 


G«  =  +  I  K-’-r*'"  -  (6.45) 

which  for  ^-*00  which  for  |t|  <<  |^|  matches  to  (4.10),  as  it  must  for  a  consistent 
matching. 
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6c.  5-2/7  ^  1 

In  this  case,  \t\  »  1.  For  |r|  >>  |^|  >>  1  ,  it  is  clear  that  the  r  dependence 
in  the  as3rniptatic  behavior  will  disappear  at  each  order  (though  not  uniformly).  Thus,  for 
large  t  the  inner  equation  (6.6)  for  r  independent  solution  becomes: 

0  =  +  ^  -  2  (6.46) 

On  substituting  '  =  H  ,  vfe  get  the  inner  equation  obtained  in  the  steady  state 

analysis  of  Combescot  et  al  (1988)  for  the  relative  finger  width  A  <  ^  .  Via  Borrel 

summation,  they  showed  that  it  is  possible  to  have  solutions  so  that 

-4  e  (6.47) 

for  Arg  ^  in  [0,  ^)  .  However,  in  the  steady  state  problem,  it  is  necessary  to  have  (6.24) 
valid  for  Arg  ^  in  [— f ,  f  ]  ,  or  equivalently  (for  symmetric  steady  finger  with  a  smooth 
tip)  require  Arg  =  0  when  Arg  ^  =  0  for  sufficiently  large  ^  .  These  were 
foimd  not  to  be  valid  and  so  a  steady  finger  was  not  possible  in  that  case.  In  our  time 
dependent  problem  with  B  independent  initial  conditions,  we  have  a  different  picture. 
Assun  e  at  each  instant  of  time,  a  zero  Co(0  move  towards  the  physical  domain 

as  fast  as  the  corresponding  daughter  singularity  ^^-ue  for  the  Sailman  family 

of  solutions.  Then  it  is  enough  to  require  that  the  solution  G(  satisfy  (6.24)  for  Arg  ^  in 
[0,  7r/2]  and  a  separate  solution  (not  an  analytical  continuation)  of  (6.46)  satisfy  (6.47) 
for  Arg  ^  in  [— 7r/2,0]  .  These  two  solutions  need  not  be  identical  since  we  can  have 
branch  cut  connecting  Co(t)  and  1  <<  ^  <<  r  <<  ,  this  corresponds  to  a  cut 

between  f  =  0  and  i  =  t  .  Thus,  for  r  »  ^  >>  1  a  matching  of  the  inner 

solution  is  possible  to  the  outer  solution  of  the  form  (4.1).  Clearly,  when  ^  >>  r  >>  1, 
the  above  analysis  becomes  invalid.  This  is  not  imexpected  as  there  is  a  singularity  of  the 
outer  perturbation  expansion  (l.ll)at^  =  r,  corresponding  to  C  =  Cai(0  moves 
away  from  the  inner  zone  around  C  =  Co(0*  then  becomes  necessary  to  introduce  an 
inner-neighborhood  around  C  =  We  do  that  in  the  next  section  for  t  =  0(1).  The 

result  remains  valid  for  <  >>  02/7  gjid  so  there  is  no  necessity  to  go  through  the  analysis 
for  r  >>  1  for  ^  —  r  =  0(1), 

Further,  when  r  >>  1,  each  of  the  -4/3  rd  singularities  Cp(0  was  born  at  Co(0) 
move  away  from  the  inner  region  around  (  =  Co(0-  This  will  be  dealt  with  in  the  following 
section  as  well. 

7.  Inner  scales  around  each  created  singularity  for  t  =  0(1) 

PVom,  the  result  (6.35)  and  the  relation  of  H,  ^  and  r  to  C  and  r  that  can  be  found 
from  (6.1)-(6.3)  and  (6.9),  it  follows  that  for  t  «  02/7^  Ao{t)  defined  in  (4.11)  scales  as 
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(Note  the  scaling  does  on  associated  with  that  particular  singularity).  We  make 
the  ansatz  that  when  C*i(0i  newly  created  singular  point  of  the  outer  perturbation 
expansion  as  noted  in  oc»,tion  4,  is  not  in  the  immediate  vicinity  of  =  1,  .Ao(t)  scales 
as  for  t  =  0(1)  as  well.  When  C*i(*)  becomes  very  close  |C|  =  1,  the  inner 

equations  themselves  become  global  integro-differential  equations  (as  will  be  seen  in  section 

8  for  special  cases)  and  in  that  case,  we  do  not  expect  the  scaling  to  survive.  On 

integration  of  (4.12),  we  get 

C,W  =  -  l\i(Ut),t)dt'  -  (7.1) 

^From  section  4,  C»i(t)  =  — ^1o(Cji(^)>  0*  Since  initially  Cp(0)  =  C*i(0)  =  Co(0)>  it  follows 
from  (7.1)  and  our  ansatz  about  the  scaling  of  Aq  that 

<p(t)  =  <..(«)  +  (7.2) 

Now,  consider  first  the  inner  region  around  C  —  C«t  (t).  It  turns  out  that  this  does 
not  include  or  overlap  with  the  inner  region  around  any  of  the  Cp(t)  Clearly,  as  pointed  out 
in  the  last  section,  for  t  »  ,  a  singular  point  C*i  of  the  outer  perturbation 

expansion  (4.1),  initiedly  coinciding  with  (^o  >  moves  towards  the  physical  domain.  We  now 
introduce  the  inner  equations  near  C*i(0  where  this  apparent  singularity  in  the  outer 
equation  is  smoothed  out.  The  appropriate  scaled  variables  in  this  case  are: 

C-U(<)  =  e~/o  9ioc(Cn(*').<')  -  (7  3J 

Ai{0)  efo  »'o((C.x (*').<')  .r)  (7.4) 


where 


r  970  (C.I  (<*),<*)  3  /*'  dti  9i0((C.i(‘i).‘i) 

-  Jo  zoT  «..«').  *■) 


Then,  to  the  leading  order,  we  find 


JHf  —  Hj/i/i/ 

The  appropriate  solution  to  Miis  that  matches  with  the  singularity  given  by  (4.8)  is: 

where  F  satisfies  the  ordinary  differential  equation: 

F  -  I  7  F'in)  =  Fm 
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It  is  clear  that  for  large  rj  ,  one  possible  behavior  is  given  by 


F  ~  (7.9) 

and  this  matches  with  (4.8).  In  order  to  ascertain  that  (7.9)  holds  over  a  big  enough  sector 
of  the  complex  plane,  we  consider  WKB  solution  to  (7.10),  which  clearly  can  be  a  linear 
combination  of  the  two  asymptotic  behavior: 

^1/2  ^  (7  11) 

We  can  ensure  the  asymptotic  behavior  (7.9)  for  Argfj,i.e.  Arg  ^  in  (— ^,^)  which 
is  a  combination  of  two  Stokes  sectors.  In  that  case,  for  Arg  fj  in  (0,  ^)  ,  including  the 
leading  order  transcendental  correction,  the  asymptotic  behavior  for  F  will  be 

F  -  ^-5/2  +  XJj  ^1/2  g<  ^  V®/*  (7  12) 

For  Arg  fj  in  (— x  >  ^) 

F  ~  ^-V2  +  £,2  ij^^^  e"‘  ^  (7.13) 


where  Di  and  D2  are  uniquely  determined  order  unity  constants.  Since  the  sector  of 
matching  is  bigger  that  tt  ,  then  according  to  suggested  matching  principle  in  section  5, 
an  inner-outer  matching  is  possible  in  this  case. 

Now  consider  the  inner-zone  aro\md  each  singularity  Cp(*)-  appropriate  to  intro¬ 
duce  inner  variables: 


C  -  c,(f)  =  X 


Hi,*)  = 


1 


(7.14) 

(7.15) 


Then  plugging  it  into  (2.17), 


the  leading  order  equation  equation  is 


where 


^  +  M(t)  Hstui  =  0 


Mr,  = 

ip  +  iiAipi*),*) 


(7.16) 


(7.17) 


In  order  to  match  with  the  outer  solution,  we  must  require  that  as  x  00  in  cer¬ 
tain  sectors  (to  be  specified)  that  correspond  to  certain  directions  towards  the  physical  ^ 
domain, 

H  Hoo  (7.18) 
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where 


(7,19) 


=  K(c,(<),()  r'''‘ 

Solution  to  the  ordinary  differential  equation  (7.16),  where  M{t)  is  merely  a  parameter, 
can  be  obtained  in  closed  form  (in  terms  of  elliptic  functions)  that  ensures  the  validity 
of  the  matching  (7.18)  for  Arg  x  in  ((2A:  —  1)  tt,  (21:  +  1)  tt)  for  any  specific  choice  of 
integer  k.  A  choice  of  integer  k  clearly  exists  so  that  inner-outer  matching  is  possible  for 
A'rg  (1  -  C/Cp(<))  in  (-’>'/2,  7r/2). 

8.  Final  stages  for  specific  type  of  singularities  of  zq 

In  the  previous  section,  otu:  analytical  evidence  suggests  a  very  rich  dynamics  that 
involve  not  only  the  alteration  of  initial  singularities  but  also  the  creation  of  new  singular 
points.  As  mentioned  earlier,  all  singularities  approach  the  physical  domain  |C|  =  1,  as  to 
the  leading  order,  their  location  evolves  according  to  the  relation  ^  =  ~9io(C(^)>  *)•  Thus 
eventually  each  of  Cj(0i  C«i(0  Cp(0  arbitrarily  close  to  |C|  =  1-  In  that  case 

the  leading  order  approximation  q\  ~=  510,  q^  =  97^,  etc.  that  were  used  in  deriving  the 
leading  order  inner  equations  in  sections  5,  6  and  7  become  invalid.  It  is  clear  then  that 
we  need  to  look  at  the  modifications  of  the  global  integral  terms  q\,  57,  etc,  created  by  the 
presence  of  singularities  that  are  very  close  to  |(|  =  1  in  order  to  study  this  late  stages. 
This  is  a  rather  involved  study  by  itself. 

On  the  other  hand,  unless  we  study  this,  our  findings  in  the  previous  sections  will 
remain  of  uncertain  relevance  to  the  physically  interesting  localized  distortions  of  the 
Hele-Shaw  interface  that  occur  in  the  presence  of  singularities  very  close  to  |C|  =  1.  In  this 
section,  we  undertake  the  task  of  deriving  appropriate  inner-equations  when  a  singularity 
^a(t)  of  ZQf  of  the  form  2ro<  ~  Bo{t)  ((  —  C«(0~^  approaches  the  physical  domain  arbitrar¬ 
ily  closely  in  the  restricted  case  ^  <  /3  <  1.  Our  main  purpose  will  be  to  show  that  z^ 
scales  as  some  inverse  power  of  B  over  this  localized  scale  implying  that  the  approach  of  the 
singularities  of  this  kind  of  the  zero  surface  tension  solution  does  indeed  correspond  to  large 
distortions  of  the  physical  interface,  which  are  smoothed  over  a  small  local  scale  by  surface 
tension  effects.  We  expect  similar  effects  by  the  approach  of  other  kinds  of  singularities 
not  considered  here  (including  the  -4/3  types  created  by  surface  tension  effects). 

In  this  case,  equation  (2.17)  is  not  a  convenient  stzirting  point,  as  many  of  the  terms 
become  of  the  same  order.  We  now  consider  (  =  e”'  to  be  a  point  on  |C1  =  1  corresponding 
to  the  physical  interface.  We  work  directly  with  (2.1)  and  (2.11)  which  remains  valid  on 
Id  =  1  provided  the  integrals  Ji  and  I4  are  taken  in  the  limit  of  (  approaching  the  circle 
from  inside.  We  introduce  the  inner-scales: 

-i(i  -  e‘vc.(())  =  er  e  (8.1) 
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(8.2) 


The  seeding  of  (6.2)  is  appropriate  in  order  to  match  with  the  outer-behavior  Z( 

(8.2) 

~  Bo(t)  (C  -  (.(*))-^ 

Further  we  introduce 

a;  =  Q(^,<) 

(8.3) 

We  introduce 

Ut)  =  R{t) 

(8.4) 

and  will  define 

p(t)  =  (R-l)  B-P^ 

(8.5) 

For  p  »  1,  we  recover  the  local  inner-equations  studied  in  section  5.  So,  we  will  assume 
p  =  0(1).  Then,  on  substituting  the  above  relation  into  (2.1),  we  find  that  to  the  leading 
order, 

(8.6) 

where  the  integral  above  is  in  the  principal  value  sense.  Now,  consider  the  kinematic 
condition  (2.11).  First,  it  is  useful  to  break  up  as 


2  f  —  ifi+oo 

J — ip— QQ 


de 


Im 


i-i'  |fr(^',t)|  [^(e,<)J  |^(e,<)l 


Im 


A(C.<)  =  /?(C,<)  +  /f(C.<) 


where 


1  /■2»  ii/  I  iv'  ■, 

/0(e“',0  =  -4  /  du'-.  ■ - — 

^  Jq  c’"  -e*"  |z^(e*v, 


t)\  7r|z((e'*'',i)| 


-  L  +  ile  u;c(e'>-',f)]  .  iie  [e**- u;c(e»•^t)] 

’  Stt/o  |^<(e‘«'',t)|"  k(e‘",<)l' 


It  is  not  difficult  to  establish  that 

- 

TTh  Jo 

-f  ~  /  du'cot 

Jo 


B^^{(  +  ip{t))  r  d  1 


TT^i 


2  |^oc(e’",t)r 

2  +  o{Bn 


2  |zoc(e‘''',t)|^ 


Again 


/«  ..  _e.-2r.+3/»,.  fi  Im  Im 

'  [«  {-«'  iff({',()i"  we.oP  . 


(8.7) 

(8.8) 

(8.9) 


(8.10) 


(8.11) 
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Using  the  scalings  (8.1),  (8.2)  and  substituting  (8.10)  and  (8.11)  into  (8.7)  using  the  relation 
for  C*/C*  derived  in  section  3  and  the  relation  (3.44),  we  find  that  it  is  necessary  to  choose 


1 

"  3(1 -;8) 


(8.12) 


to  get  a  leading  order  inner-equation  that  includes  the  curvature  and  which  can  match  to 
the  outer  solution.  Obviously,  this  scaling  is  invalid  for  ^  >  1,  since  the  assumed  small 
inner-scale  is  then  inconsistent  with  (8.1).  Further,  the  approximation  (8.10)  become 
invalid  for  <  |.  With  the  choice  in  (8.12),  we  get  (2.11)  (after  differentiating  with 
respect  to  C)  to  reduce  to 


where 


H,  =  i  £-'“•«)  H 

Jo 


2  koc(e‘",<)l' 

v'  —  Va(t)  d  1 


du'  |^oc(e*'",f)| 


1 

=  +37  / 

ir'i  Jd 

A  -  fi  Ji!-  X 

'  '  U"  y-if-oo 


(8.13) 

(8.14) 
(8.16) 

(8.16) 


It  is  then  necessary  to  solve  the  system  of  integro  differential  equations  (8.6)  and  (8.13), 
subject  to  the  asymptotic  behavior  that 


~  Bo(<)HC.(<)]"^  (8.17) 

and  for  early  enough  times,  when  p{t)  »  1,  we  want 

~  5o(<)l-iC.(0]“^  (8.18) 


for  edl  We  have  not  solved  these  equations.  Nonetheless,  the  scalings  obtained  in 
(8.12),  together  with  the  relations  (8.1)  and  (8.2)  immediately  imply  that  a  large  localized 
distortion  occurs  in  the  physical  interface  at  a  point  corresponding  to  ^ 

Thus,  we  can  conclude  that  an  initial  singularity  of  Z(^(C,  0)  atleast  of  the  type  (C  —  C«(0))~^ 
with  ^  in  the  range  (^,1)  initially  located  at  0(1)  distance  from  =  1,  while  being 
smeared  out  and  transformed  by  surface  tension  effects,  will  still  make  its  presence  felt  on 
Id  =  1  after  0(1)  time  that  can  be  computed  by  calculating  the  trajectory  of  d(f)  that 
moves  according  the  the  zero  surface  tension  solution  and  eventually  comes  arbitrarily 
close  to  Id  =  1. 
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Conclusion 

We  have  imbedded  the  zero  surface  tension  problem  in  a  well  posed  problem  and  shown 
how  surface  tension  affects  singularities  of  certain  forms.  Here  we  only  addressed  the  global 
qualitative  features  of  the  dynamics  of  singularities  emd  zeros.  These  will  be  supplemented 
with  actual  numerical  calculations  of  singularity  trajectories  relating  it  to  the  shape  of  the 
physical  interface.  We  also  need  to  account  for  the  final  stages  of  singularities  other  than 
the  special  types  considered  in  section  8. 

However,  in  order  to  completely  relate  the  findings  with  physical  domain  initial  con¬ 
ditions,  which  in  general  cannot  be  analytically  continued  exactly  to  the  complex  plane, 
one  needs  to  address  the  issue  of  the  probability  density  of  singularities  of  different  types 
and  their  strengths  in  |C|  >  1,  once  a  function  is  specified  in  |C|  <  1  within  a  certain  error 
(say  in  the  Max  norm  sense).  This  does  not  make  sense  at  this  stage  since  the  singuleirity 
types  are  not  even  denumerable.  However,  if  the  long  time  feature  of  the  interface  is  not 
crucially  dependent  on  the  precise  singularity  nature,  it  may  meike  sense  to  restrict  to  spe¬ 
cific  kinds  of  singularities.  Our  findings  in  this  paper  is  in  some  sense  equivalent  to  finding 
the  input-output  relationship  of  a  device.  However,  in  experiment  or  in  direct  numerical 
simulation  (if  one  is  possible  for  very  small  H)  what  one  observes  is  the  statistical  features 
of  the  output,  given  random  inputs,  the  randomness  coming  in  because  one  cannot  control 
the  location  of  singularities  in  the  unphysical  domain  even  with  arbitrarily  accurate  data 
in  the  physical  domain.  Thus  to  predict  the  statistical  features  of  the  output,  besides 
figuring  out  the  input-output  relationship  that  has  been  addressed  to  in  this  paper,  we 
also  need  to  predict  the  statistical  features  of  the  input.  We  cein  also  expect  that  such 
statistical  features  of  the  input,  combined  with  the  understanding  of  the  dynamics  with 
specific  input  is  likely  to  lead  us  to  discover  self-similarities  of  the  statistical  features  in 
time  that  may  well  explain  the  apparently  fractal  nature  of  the  interface  over  some  range 
of  length  scales. 

We  can  perhaps  expect  this  mode  of  attack  to  be  useful  for  the  understanding  of  the 
time  dependent  features  (such  as  side-branching)  of  a  needle  crystal,  since  the  importemt 
element  of  the  analytical  technique  is  not  the  conformal  map,  but  but  a  system  of  integro 
differential  equations,  which  is  available  in  that  case. 
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y  =  -a 


2.  Fig.  2:  Rectilinear  Hele-Shaw  flow  viewed  in  the  lateral  plane. 
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4.  Fig.  4:  Hele-Shaw  flow  in  a  radial  geometry. 
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6.  Fig.  6:  Illustration  of  matching  zone.  TranscendenteJly  large  terms  in  B  have  to  be 
avoided  inside  circle  of  radius  R  where  the  nearest  singularity  S  is  just  beyond  R  . 
The  sector  of  matching  near  5  is  bounded  by  the  two  dotted  lines  through  S  and 
subtend  em  angle  of  . 


9.  Fig.  9:  Branch  cut  between  ^  and  corresponding  daughter  singularity 

^  allows  two  different  local  solutions  on  two  sides  of  the  cut  to  be  matched  to 

the  outer  solution  in  the  region  inside  the  large  dotted  circle. 
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